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FOREWORD 

The  work  reported  herein,  covering  the  period  10  April  1972  to  31  December 
1972,  was  carried  out  by  the  Infrared  and  Optics  Division  of  the  Environmental 
Research  Institute  of  Michigan  (formerly  the  Willow  Run  Laboratories  of  The 
University  of  Michigan),  Ann  Arbor,  Michigan.  The  work,  which  was  performed 
under  Contract  DAAD05-72-C-0246  for  the  Army  Ballistic  Research  Laboratories, 
was  done  in  three  parts,  each  of  which  represents  one  volume. 

The  three  volumes  comprising  the  present  series  on  polarized  radiance  are: 

I  -  Polarized  Bidirectional  Reflectance  with  Lambertian  or  Non-Lambertian 
Diffuse  Components 

II  -  Polarized  Spectral  Emittance  from  4  to  14  nm 
III  -  Wavelength  Dependence  of  Polarized  Bidirectional  Reflectance 
The  ERIM  number  for  Volume  I  of  this  report  is  192500- l-T(I). 


PREFACE 


Under  contract  commitments  with  the  Air  Force  Avionics  Laboratory  (Con¬ 
tract  F336 15 -70-0 1123),  a  bidirectional  reflectance  model  was  developed  to 
predict  reflectance  from  rough  surfaces  as  a  function  of  the  zenith  and  azimuth 
angles  of  source  and  receiver,  given  a  set  of  input  parameters  derived  from 
limited  measurements.  It  was  observed,  however,  that  the  model  required 
further  extension  to  account  for  anomalies  in  the  comparison  of  measured  data 
with  calculated  model  predictions.  Specifically,  it  was  determined  that  the  model 
would  be  more  accurate  if  it  could  account  for  a  non-Lambertian,  non-specular 
component  of  the  reflectance  which  was  assumed  to  result  from  scattering  within 
the  target  material.  Extension  of  the  model  to  take  the  non-Lambertian  angular 
dependence  into  account  was  carried  out  under  this  BRL  contract,  as  well  as 
coding  the  model  in  Fortran  IV  and  validating  it. 

This  extension  of  the  bidirectional  reflectance  model  plus  other  extensions 
performed  under  the  AFAL  contract  have  considerably  improved  the  fit  between 
model  predictions  and  measured  data. 

The  work  done  under  this  contract  with  BRL  has  been  combined  with  the 
work  done  under  the  Air  Force  contract  to  form  a  unified  model.  Therefore, 
for  completeness,  the  AF-sponsored  part  of  the  model  is  included  in  this  re¬ 
port. 
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POLARIZED  RADIANCE 
Volume  I 


Polrrized  Bidirectional  Reflectance  with  Lambertian 
or  Non-Lambertian  Diffuse  Components 

10  April  Through  31  December  1972 
1 

INTRODUCTION 

A  model  for  predicting  the  radiance  at  a  remote  sensor  must  include  the  spatial,  spectral, 
and  polarization  characteristics  of  the  bidirectional  reflectance  and  directional  emittanee  with 
respect  to  target  and  background  surfaces.  In  principle,  the  directional  reflectance  and  di¬ 
rectional  emittanee  properties  of  materials  must  be  known  for  all  source  and  receiver  angles, 
polarizations,  and  wavelengths.  A  Lambertian  assumption  may  be  valid  for  some  types  of 
backgrounds,  but  for  most  man-made  targets  is  scarcely  adequate.  Measurement  of  all  spatial, 
polarization,  and  spectral  characteristics  of  the  bidirectional  reflectance  and  directional  emit- 
tance  for  a  large  number  of  material  samples  is  impractical.  Even  if  such  measurement  were 
performed,  the  data  could  not  all  be  stored  efficiently  enough  to  make  it  accessible  for  digital 
computations.  Clearly,  an  empirical  model  is  required  to  approximate  the  bidirectional  reflec¬ 
tance  and  directional  emittanee  properties  from  a  limited  number  of  measurements. 

The  bidirectional  reflectance  model  first  developed  by  the  Environmental  Research  Insti¬ 
tute  of  Michigan  (ERIM)  for  the  Air  Force  [1,  2]  is  described  in  this  report.  The  model  accounts 
for  effects  that  produce  both  specular  and  diffuse  components.  In  particular,  a  surface  model 
relates  bidirectional  reflectance  for  all  source -receiver  angles  and  polarizations  to  fixed 
bistatic  measurements  and  a  Brewster  angle  measurement.  The  model  has  been  extended 
under  this  contract  to  enable  calculation  of  either  a  Lambertian  diffuse  component  or  a  non- 
Lambertian  diffuse  component.  The  latter  component  accounts  for  angular  and  depolarization 
properties  arising  from  internal  scattering  effects.  Our  extension  of  the  bidirectional  re¬ 
flectance  model  has  considerably  improved  the  fit  between  model  predictions  and  measured  data, 
as  will  be  shown  in  Section  6. 

As  it  now  stands, the  model  permits  generation  of  an  enormous  amount  of  bidirectional  re¬ 
flectance  data  from  a  very  small  amount  of  measured  data.  The  accuracy  shown  in  Section  6 
on  Model  Validation  indicates  that  the  model  is  very  effective,  although  it  can  still  be  improved, 
particularly  at  large  receiver  zenith  angles.  With  the  ability  to  account  for  elliptical  (partic¬ 
ularly  circular)  polarization  now  built  in,  the  model  is  available  for  use  with  circularly  po¬ 
larized  sources,  if  these  sources  prove  useful  in  the  future. 

In  this  report,  we  compare  measured  data  with  results  computed  from  both  the  initial 
model  and  from  the  extended  model, and  then  evaluate  the  relative  performance  of  the  two.  We 
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establish  a  domain  of  validity  for  each,  based  on  material  properties.  Since  the  modeling  is 
empirical,  only  a  limited  amount  oi  measured  data  are  required  as  input  parameters.  In  this 
case,  the  parameters  are  the  fixed  bistatic  data. 

All  modeling  described  in  Volume  I  of  this  report  was  performed  with  respect  to  one  wave¬ 
length,  A  =  1.06  4 m. 
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2 

BIDIRECTIONAL  REFLECTANCE 


One  physical  property  which  can  be  measured  directly  from  a  sample  of  material  is  bi¬ 
directional  reflectance.  The  physical  definition  is 


cLr(tfr,  0r) 

6E.(0.,  0.) 


(la) 


where  <5E.(0.,  0.)  is  the  incremental  irradiance  (power  per  unit  area)  impingent  on  the  surface 
of  a  material  fromthe direction  0.),  and  Lr(f?r>  </>r)  is  the  resulting  increment  of  radiance 
(power  per  unit  projected  area  per  unit  solid  angle)  scattered  from  that  surface  in  the  direc¬ 
tion  (<9r,  0r).  Figure  1  illustrates  the  situation.  The  bistatic  angle,  2)3,  is  that  angle  between 
the  vectors  which  point  to  the  source  and  the  receiver  respectively. 

Equation  (la)  can  be  rewritten  in  terms  of  directly -accessible  experimental  parameters  as 

6P 

_ r 

i>A  cos  9  6i1 

0  )  = - jp1 - L  (lb) 

_ i 

OA 


where  t>P.  is  the  power,  inwatts,  incident  from  the  direction  (9.,  0.)  on  the  small  area  <5A,  and 

6P  is  the  resulting  power  scattered  into  the  small  solid  angle  69  in  the  direction  (9  ,  <t>  ). 
r  t  vr  r 

When  polarization  dependence  is  to  be  shown,  subscripts  are  appended  tothep'  term.  Thus 

when  we  write  p’  ,  the  leading  subscript,  a .,  describes  the  source  polarization  while  the  trail¬ 
er  1 

ing  subscript,  £*r,  describes  the  receiver  polarization.  The  source  polarization,  always  re¬ 
ferred  to  the  plane  of  incidence,  describes  the  polarization  state  of  the  electric  field  vector. 

The  appended  subscript  symbols  II  and  1  indicate  whether  the  source  electric  vector  polariza¬ 
tion  is  parallel  to  or  perpendicular  to  the  incidence  plane.  The  reflected  electric  field  polariza¬ 
tion  state  is  specified  by  the  same  symbols,  but  here  the  reference  plane  is  that  reflectance 
plane  defined  by  the  sample  normal  and  the  direction  to  the  receiver.  (For  example,  p!  ^  re¬ 
presents  reflectance  measured  when  source  polarization  is  perpendicular  to  the  incidence 
plane  and  receiver  polarization  is  parallel  to  the  reflectance  plane.)  Notice  that  when  either 
the  source  or  the  receiver,  or  both,  are  scanned  in  angle  over  the  sample,  the  incidence  and 
reflectance  planes  change  orientation  with  relation  both  to  the  sample  and  to  each  other. 

Bidirecticnal  reflectance  depends  on  the  physical  properties  of  the  material  as  well  as  on 
the  geometric  state  of  its  surface.  Different  surface  states  result  in  different  reflectances. 
Hence,  a  complete  collection  of  bidirectional  reflectance  data  for  any  single  material  would 
require  measurements  of  a  large  number  of  samples  of  the  material, each  with  a  different  sur¬ 
face  state.  Each  sample  would  have  to  be  measured  with  several  source -receiver  polarization 
combinations.  Consequently  a  very  large  number  of  source  and  receiver  positions  would  be 
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required  for  each  set  of  polarization  states.  Finally,  the  entire  procedure  would  have  to  be 
repeated  at  many  different  wavelengths.  The  purpose  of  modeling  is  to  predict  reflectance 
data  from  only  a  limited  number  of  measurements  and  hence  eliminate  the  need  for  an  other¬ 
wise  unwieldy  measurement  program. 
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BACKGROUND  INFORMATION 


The  bidirectional  reflectance  model  is  based  on  observations  of  polarized  bidirectional 
data  from  rough,  painted  surfaces  which  exhibit  a  Brewster  angle  (Fresnel -like  behavior  in 
relation  to  the  specular  geometry).  The  degree  of  depolarization  appears  slight,  and  on  that 
basis  for  the  initial  modeling  work  in  this  program  single  specular  reflection  from  the  rough 
front  surface  was  assumed  to  be  the  dominant  reflection  mechanism.  Multiple  front -surface 
reflections  and  internal  scatterings  were  observed  to  be  smaller  and  were  initially  incorpo¬ 
rated  into  a  Lambertian  "volume"  model  to  account  for  the  diffuse  component. 

The  assumption  that  the  diffuse  component  is  Lambertian,  however,  makes  it  difficult  to 
account  for  certain  anomalies  that  occur  when  measured  data  are  compared  with  the  model’s 
output.  For  example,  Fig.  2  is  a  bidirectional  reflectance  curve  showing  the  reflectances  at 
the  receiver  as  the  receiver  scans  over  zenith  angles  from  0°  to  90°  in  the  0^  =  0°  and  <t>T  =  180° 
half  planes.  The  source  remains  fixed  at  y.  =  40°  and  0.  =  180°.  The  upper  curve  shows  re¬ 
flectances  when  source  and  receiver  are  both  linearly  polarized  at  the  same  polarization  angle 
with  respect  to  the  target -incidence  and  target -receiver  planes.  (In  this  case,  both  are  per¬ 
pendicular-polarized.)  The  lower  curve  shows  reflectances  when  source  and  receiver  are 
cross -polarized  with  respe  t  to  one  another.  (Source  is  perpendicular -polarized;  receiver  is 
parallel -polarized.)  Note  the  marked  angular  dependence  in  the  lower  curve.  If  the  nonspecular 
component  were  truly  Lambertian,  no  such  angular  dependence  would  be  present. 

Also,  although  radiation  sources  in  this  work  are  all  linearly  polarized,  future  work  may 
well  involve  more  general  cases.  Therefore,  the  model  should  account  for  the  most  general 
type  of  polarization  — namely,  elliptical. 

For  the  above  reasons,  and  in  order  to  obtain  a  closer  overall  correspondence  between 
model  prediction  and  measured  data,  the  model  has  been  extended  to  account  for  the  following: 

(1)  possible  non-Lambertian  angular  dependence  of  depolarized  component 

(2)  shadowing  and  obscuration  produced  by  the  roughness  of  the  surface 

(3)  elliptical  polarization 

The  model — a  phenomenological  one  in  that  its  use  requires  a  limited  number  of  measure¬ 
ments  — is  described  in  the  next  two  sections.  Section  4  includes  a  discussion  of  specular  re¬ 
flectance  from  the  surface,  effects  caused  by  shadowing  and  obscuration  resulting  from  surface 
roughness,  and  polarization  effects.  Section  5  describes  the  volume  model. 
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BIDIRECTIONAL  REFLECTANCE  (STER.'1) 
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SURFACE  MODEL 

In  this  section,  we  review  the  surface  model  and  also  discuss  the  interference  effects  that 
necessitated  model  modification. 


4.1.  AVAILABLE  AREA 

If  the  rough  surface  is  considered  to  be  made  up  of  small  sequins  having  a  distribution  of 
orientations,  there  will  be  some  specular  reflectance  at  any  receiver  angle  and  the  extent  of 
that  reflectance  will  be  determined  in  part  by  the  amount  of  surface  oriented  for  specular  re¬ 
flection  at  that  receiver  angle.  (The  area  available  for  such  reflection  will  also  depend  on  how 
some  sequins  shadow  or  obscure  others.)  Since  measurements  do,  in  fact,  show  a  reflectance 
distribution  over  the  hemisphere,  we  assume  that  the  above  description  is  valid  and  that  there  is, 
indeed,  a  distribution  of  surface  areas  which  have  normals  pointing  in  different  directions.  There¬ 
fore,  to  establis.  •  ■  -  iis'ribution  of  available  surface  area,  we  define  a  density  function  H(0,  <t>) 
which  describes  the  re' Jive  density  of  local  surface  normals  (per  steradian)  pointing  in  the 
direction  (6,  <t> ). 

The  effect  of  the  distribution  of  surface  normals  is  measured  by  a  zero  bistatic  measure¬ 
ment  in  which  <9.  =  «  and  *•.  -  <t>„-  (Note  that  we  really  use  a  fixed  bistatic  scan  with  a  small 
l  r  l  r 

bistatic  angle.  A  true  zero  bistatic  scan  would  be  very  difficult  to  obtain  since  source  and 
receiver  obviously  cannot  occupy  the  same  position.) 


4.2.  FRESNEL  COEFFICIENTS 

Fresnel  reflectance  coefficients  describe  the  reflectance  and  polarization  of  specularly 
reflected  radiation  as  functions  of  source  and  detector  positions  and  of  the  complex  index  of 
refraction.  However,  since  we  are  trying  to  find  reflectance  as  a  function  of  source  and  de¬ 
tector  positions  only,  we  must  know  — or  be  able  to  determine  — the  index  of  refraction.  (As 
discussed  later  in  this  section,  we  can  determine  the  index  by  measuring  the  Brewster  angle.) 
Since,  in  the  surface  model,  we  consider  only  single,  local  specular  reflections,  the  Fresnel 
equations  automatically  account  for  polarization. 


If  the  receiver  subtends  the  solid  angle  6flr  from  the  sample  (see  Fig.  1)  the  solid  angle 
t>nA  in  which  local  surface  normals  must  lie  to  permit  collection  of  the  local  specularly  re¬ 
flected  radiation  by  the  receiver  is  given  by: 


6«r 

—  cos  (i 

ft  4 

This  solid  angle  is  centered  about  the  direction  (a  .  ^  V 

'  n  ft' 


(2) 
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Let  £>P.  be  power  incident  on  area  <>A.  The  fraction  of  surface  area,  6A[y  ,  ^  which 


reflects  radiation  into  the  receiver  is  given  by 


ft  ft 


W0 


0  =  E  ,  <t>  <3A 


A  A  i 

n  n 


ft 

The  power  incident  on  <5A 


ft 


(v 


n  n 


(3) 


OP 


6A(fl  ,  <t> 

ft  cos  /3 


<5A 


cos  9. 

l 


(4) 


Since  the  Fresnel  reflectance,  R(/3),  is  just  the  ratio  of  reflected  power  to  incident  power,  then 

6A(e  ,  <t>  \ 

dP^R03)6P, 


cos  9 


(5) 


<3  P_ 


6  A  cos  6  6S1 

Recall  that  in  Eq.  (lb):  p'(9.,<t>.\  9  ,<t>  )  =  - ttet - “•  Substituting  Eqs.  (5),  (3)  and  (2): 

i  l  r  r °  v 


9  <f>  )  = 


R (0)H(6  .  0  ] 
V  ft  ft/ 


<5P. 

“OA 


i,9i’  r,?V  4  cos  cos  0 
i  r 


(6) 


By  considering  the  case  when  source  and  receiver  are  in  the  same  position,  i.e.,  a  zero 
bistatic  (fl  -  0 )  case,  E(  9  ,  0  ]  can  be  determined.  In  this  situation 

V  ft  ft/ 


R(0)  E  ( 9  ,  0 


p'(o  .0  ;  a  .«  )  = - ^ 

'  ft  ft  ft  n '  4  cos 


ft 


and 


ft 


2(4,  0) 


4 p'(o  ,<t>  ;  9  ,<t>  ]  cos^  (9 

'  ft  ft  ft  ft/  ft 


"W 


(7) 


(8) 


Now  substitute  back  into  Eq.  (6)  and 

p'(  a  ,<t>  ;  9  ,<t>  'j 
RCi)  '  ft  ft  ft  ft  / 


2 

cos  9 


f)  (  .  0 .  J  9  ,  0  )  -  yr 

*  '  i  i’  r  r  R(0 


cos  9.  cos  9 
l  r 


(9) 


Equation  (9)  is  an  expression  for  the  bidirectional  reflectance  given  in  terms  of  measured 
data  and  Fresnel  reflectance  coefficients.  However,  to  evaluate  the  Fresnel  coefficients  so 
they  can  be  used  in  Eq.  (9)  takes  a  little  work.  For  example,  R(/I)  is  a  function  of  the  real  and 
Imaginary  parts  of  the  complex  index  of  refraction,  n’  =  n  -  ik  (see  Ref.  3  or  Appendix  in). 
Therefore,  n  and  k  must  be  found  before  R(/3)  can  be  determined. 
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Moreover,  k  is  taken  to  be  very  small*  so  that  n  can  be  determined  experimentally  by  mea¬ 
suring  the  Brewster  angle,  t)  and  then  using  n  tan  f/_  to  solve  for  n. 

o  B 


4.3.  SHADOWING  AND  OBSCURATION 

Referring  to  Eq.  (9),  we  can  derive  a  zero  bistatic  curve,  p’(PA  ,4>  ;  0A,tf>A  from  a  p'{9  <t> 

\  n  ri  n  n/  i  i 

<>  0  )  curve  with  9  </>.  lixed  and  w  variable  by  inverting  the  equation  so  that 


9  4>  ,  9  <p 

A  A  A  A 

n  n  n  n 


R(0) 

W) 


p'(f>.0.  9  <t>  )  cos  9.  cos  9 


l  l  r  r 


2 

COS  I) 


(10) 


A 

n 


after  doing  this  for  a  variety  of  fi.’s,  we  found  that  the  curves  obtained  differed  systematically 
from  those  obtained  from  a  fixed  bistatic  measurement.  Apparently,  because  of  surface  rough¬ 
ness,  some  sequins  shadow  or  obscure  others;  this  reduces  reflectance  everywhere  except  at 
a  purely  back -scattered  position.  The  model  must  therefore  be  modified  to  correct  for  such 
interference. 

Torrance  and  Sparrow  [4]  have  developed  an  analytical  function  that  helps  correct  the  sit¬ 
uation;  however,  we  have  constructed  our  own  function  using  empirical  considerations  only. 

Our  function  results  in  better  agreement  between  measured  and  derived  fixed  bistatic  curves 
than  does  the  analytical  function  of  Torrance  and  Sparrow.  The  empirical  function  (SO)  is  de¬ 
fined  as: 


SO  - 


e . 

_n  -2/?/t 
SI  l 


6, 


1  * 


%  A 


<t>  0. 

_n  _i 


(ii) 


where  SI  and  r  are  parameters,  and  <t>  is  a  factor  calculated  from  the  geometry,  which  adjusts 
the  fall-off  rate  of  the  shadowing  and  obscuration  function  in  the  forward -scattered  direction. 


*For  the  calculations  in  this  study,  results  of  past  measurement  programs  [1]  were 
used  to  establish  the  refractive  indices.  In  those  programs,  it  was  determined  that  the 
magnitude  of  the  total  index  of  refraction  was  close  to  1.65;  that  the  imaginary  part  of  the  in¬ 
dex  of  refraction  could  be  neglected,  compared  to  the  real  part;  and  that  the  index  of  refraction, 
for  the  wavelengths  of  incident  radiation  under  consideration  (1  to  4  pm),  did  not  vary  apprecia¬ 
bly. 
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We  now  modify  Eq.  (9): 


P'O^J  or,<t>r)  - 


p'(l>  ,<>  ;  h  ,<p  )  cos2  " 

R(i3)  '  ft  ft  ft  ft ' _ ri 

R(0)  cos  t).  cos  i) 

1  r 


(SO) 


Equation  (10)  becomes 


,<t>  ;  6  ,<f> 

ft  ft  ft 


yr.^r)  cos  6.  cos  Pr 
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VOLUME  MODELS 

The  following  discussion  outlines  the  reasoning  behind  the  extended  portion  of  the  bidirec¬ 
tional  reflectance  model.  (The  extended  portion  is  referred  to  as  the  "volume"  model.) 

Different  materials  with  varying  degrees  of  surface  roughness  and  different  optical  prop¬ 
erties  show  differences  in  nonspecular  reflectance  behavior.  These  differences  show  up  in  the 
extent  to  which  the  nonspecular  reflectance  is  dependent  upon  angular  position  of  the  receiver. 

To  make  provision  for  materials  that  do  exhibit  such  angular  dependence  and  for  those  that 
do  not,  two  volume  models  are  used.  The  following  discussion  describes,  first,  a  Lambertian 
volume  model  which  has  no  angular  dependence,  and  then  a  non -Lambertian  volume  model  in 
which  angular  dependence  is  important. 

5.1.  LAMBERTIAN 

In  addition  to  Fresnel  reflection  from  a  surface,  other  effects  such  as  might  take  place  be¬ 
neath  the  surface  can  produce  a  nonspecular  reflectance  component  everywhere  in  the  hemi¬ 
sphere.  If  the  surface  roughness  as  well  as  the  absorption  properties  of  the  surface  are  right, 
this  volume  reflectance  may  be  completely  diffuse  and  uniform  over  the  hemisphere.  Moreover, 
the  reflected  radiation  will  be  totally  depolarized,  regardless  of  the  polarization  of  the  source. 
Thus,  if  the  receiver  is  polarized  in  the  orthogonal  direction  to  the  source  polarization,  an  in¬ 
plane  measurement  will  represent  the  volume  component  only.  However,  only  half  the  volume 
component  is  actually  represented,  since  there  should  be  an  equal  diffuse  contribution  polarized 
in  the  same  direction  as  the  source. 

The  Lambertian  volume  component  is  one  of  the  input  parameters  for  the  model  when  a 
target  material  with  Lambertian  reflectance  properties  is  considered.  A  method  whereby  values 
for  this  parameter  may  be  extracted  is  described  in  Section  6. 

5.2.  NON -LAMBERTIAN 

On  the  basis  of  the  Lambertian  diffuse  model  described  above,  no  angular  dependence 
would  be  expected  for  the  diffuse  component.  However,  for  some  materials,  actual  measure¬ 
ments  show  that  there  is  an  angular  dependence.  To  provide  for  the  angular  dependence  of  the 
diffuse  component,  the  model  has  been  extended  by  including  scattering  that  takes  place  beneath 
the  surface. 

Assuming  an  exponential  scattering  function  as  the  radiation  first  enters  and  then  leaves 
the  surface,  and  making  reference  to  Fig.  3,  we  construct  an  expression  for  the  volume  scat¬ 
tering  component  of  the  bidirectional  reflectance  as  follows: 
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Ei0  irradiance  at  surface  of  area  A^,  where  A^  is  the  area  of  cross  section  of  the  colli¬ 
mated  beam  and  is  normal  to  the  beam 
Eq  =  irradiance  on  surface  element  of  area  A 

E  irradiance  on  surface  ol  slab  or  area  A  at  distance  z  beneath  surface 
i 

-  radiance  scattered  from  primary  beam  through  20  in  direction  of  receiver 
ji  ^  half  of  angle  between  target -to -source  and  target -to-receiver  vectors 
a  =  total  scattering  cross  section  (ignoring  absorption) 
o(d)=  differential  scattering  cross  section  with  respect  to  (1,  i.e.,  / o  (/ 3 ) df2  =  /da/dfldtt  =o 
fl  ^  solid  angle  subtended  at  target  by  receiver,  assuming  a  point  target 
•/.  -  angle  of  incident  beam  relative  to  fixed  z  axis 
=  angle  of  reflected  beam  relative  to  fixed  z  axis 

The  objective  of  the  following  calculation  is  to  determine  that  portion  of  the  primary  beam 
scattered  from  distance  z  beneath  the  surface  through  an  angle  2 0  toward  a  receiver  which  sub¬ 
tends  solid  angle  $2 . 


First,  the  bidirectional  reflectance  defined  in  Eq.  (1)  is  now  p'  =  L r/Eg  with  respect  to  the 
slab  (see  Fig.  3).  To  determine  E^: 


A  i  -  A  cos  a. 


E  .  E _ _P _ 

0  A  A,/ cos  n. 


P  cos  «. 
l 


Ho 


(14) 

(15) 


where  P  is  the  power  at  surface  of  area  A. 

The  irradiance  incident  on  the  slab  at  distance  z  beneath  its  surface  is: 
t  -az/cos  a. 

Ej  EqO  E0e  (16) 

where  (  z  cos  Hence 

i 

-a  z  cos  a, 

E.  --  Ej0  cos  Ojt*  (17) 

and 

-ai  cos 

dEj  dz  (18) 

where  dE,  is  the  amount  by  which  irradiance  decreases  in  going  from  distance  z  to  distance 

-o  z  cos  -oz  cos  n 

z  *  dz  beneath  surface.  Note  also  that  e  1  and  e  r  represent  the  scattering 

loss  from  the  beam  on  the  way  in  and  on  the  way  out  of  the  material,  respectively.  To  deter¬ 
mine  L  . 
r 
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Pr  =  LrAJlr  cos  0r  -  power  at  the  receiver  (19) 

dLf  =  radiance  scattered,  in  direction  of  receiver,  from  one  small  slab  (20) 

of  thickness  dz 

Radiance  from  slab  in  direction  9  (or  fl)  can  be  written: 

r 

dLr  =  -dE.o(0)  (21) 

since  o(/3)  is,  by  definition  the  fraction  of  beam  scattered  into  2/3 .  Note  that  (since  we  are  ignor¬ 

ing  absorption)  irradiance  lost  from  the  incident  beam  is  the  radiance  of  the  scattered  beam; 
therefore  a  minus  sign  precedes  dE..  Hence,  if  there  is  no  further  power  loss 

dPr  =  -A£lrdE.cr(i3)  (22) 


However,  power  loss  caused  by  beam  scattering  occurs  on  the  way  out  as  well  as  the  way  in; 

-oz/cos#^ 

on  the  way  out. 


the  loss  is  represented  by  e 


Therefore 


dP 


-oz/cos  9 


=  -A«  o(/i)e 
r  r 


dE.  =  -- 


i  /  cos  0.  r 


n  a(j3)e 


-oz/cos 9 


dE. 


(23) 


Substituting  the  expression  for  dE.,  Eq.  (18),  into  Eq.  (23),  we  obtain: 

/  -0Z/COS0 

A  ,o| 


-oz/cos  0, 


dP  =  E,  noe 
r  10 


r(/3)e 


cos  9, 


i 


)  «r  das 


Pr  r 


E10AlS1r 
dP  -  — L  £ 
r  cos  0 


o(/j) 

\cos  0j  cos  9J 


(24) 


(25) 


where  the  Integration  from  0  to  <x>  assumes  no  transmission  of  power  through  the  material,  l.e., 
the  material  has  effectively  an  infinite  thickness  with  respect  to  transmission.  Therefore 


Ei0o(/3) 


r  A  cos  9JI  77777  I  I  V 

r  r  cos  9  - —  +  - — 

r\cos  Oj  cos  0  / 


and 


P‘  = 


"0  cos  9,  cos  9 


_  o(d) 

1  \  "  (cos  9.  +  cos  9  ) 


l  r\,cos  0^  cos  9 


3 


i 


(26) 


(27) 


In  ignoring  the  finite  thickness  of  the  layer  of  material,  we  have  also  ignored  the  possible 
specular  reflectance  of  the  bottom  surface.  To  account  for  the  possibility  of  specular  reflec¬ 
tion  from  the  bottom  layer,  it  may  be  useful  to  provide  a  parameter  function  peaked  near 


A  =  0.  Therefore,  we  include  all  ji  dependence  in  a  function  f(/3),  and  all  0A  dependence  in  a 
n 


function  g 


k).* 

'  n ' 


'A 

n 


and  write 


P'  =  2 


Pyf(0)g^  J 

cos  0.  +  cos  9 
i  r 


(28) 


where  f(/l)  and  g((9 A  j  provide  freedom  for  empirical  adjustment.  The  constant,  p^.,  rep-esents 


n  /  \ 
the  value  of  p'  when  «.  =  =  0  and  f{/3)  =  g(flA  j  =  1. 
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MODEL  VALIDATION 

Use  of  the  bidirectional  reflectance  model  requires  a  limited  amount  of  measured  data 
(namely  the  zero  bistatic  measurement)  from  which  complete  sets  of  reflectances  can  be  cal¬ 
culated.  The  results  of  these  model -calculated  bidirectional  reflectances  can  then  be  com¬ 
pared  to  corresponding  results  of  actual  measurements.  This  was  the  procedure  we  followed 
to  validate  the  model. 

Model  calculations  and  measured  data  were  compared  in  terms  of  p'  (the  reflectance),  a  or 
i/'  (the  angle  of  polarization  for  the  beam  after  reflection  from  the  target),  and  P  (the  percent¬ 
age  of  polarization  of  the  reflected  beam). 

Measured  data  for  materials  of  different  properties  (color  and  roughness)  were  used  to 
demonstrate  the  model's  performance.  The  materials  are  designated  as  A02016-001, 
A02018-002,  and  A02100.  Material  A02018-001  is  a  green  paint  and  mater’al  A02018-002  is  a 
tan  paint.  These  materials  were  supplied  by  the  Army  Ballistic  Research  Laboratories  for  the 
purpose  of  developing  the  non -Lambertian  diffuse  component  of  the  model. 

Measured  data  for  materials  A02018-001  and  A02018-092  were  used  for  the  fitting  since 
it  was  felt  that  two  surfaces  of  extremely  different  properties  (color  and  roughness)  would  be 
necessary  to  demonstrate  the  performance  of  the  model.  Measurements  on  material  A02017- 
001  show  that  the  bidirectional  reflectance  very  closely  resembles  that  of  material  A02018-001. 
Therefore,  sufficient  information  was  developed  in  validating  the  model  with  material  A02018- 
001  to  permit  assignment  of  parameters  to  material  A02017-001  as  well.  Additional  validation 
was  performed  with  respect  to  A02100  (soil)  and  discussion  is  included.  Model  parameters  are 
listed  in  Table  I.  (See  Section  7  for  definitions  of  model  parameters.)  The  overall  discussion 
of  the  model  fitting  is  divided  into  five  parts: 

(1)  p' for  A02018-001 

(2)  p'  for  A02018-002 

(3)  polarization  angle  (a  or  \l/  )  for  A02018-001 

(4)  percent  polarization  (P)  for  A02018-001  and  A02018-002 

(5)  p'  for  A02100 

In  what  follows,  the  orientation  of  the  source  polarizer  in  the  measurements  of  materials 
A02018-001  and  A02018-002  was  not  actually  perpendicular,  parallel,  nor  at  45°  to  the  plane  of 
incidence  but  instead  was  offset  by  5°  in  each  case.  Specifically,  the  appropriate  correspon¬ 
dences,  shown  in  Table  II,  should  be  recognized.  These  shifts  were  taken  into  account  when 
the  validation  calculations  were  made  on  the  computer;  however,  we  continue  to  refer  to  "per¬ 
pendicular,"  "parallel,"  and  "45°. " 
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TABLE  I.  MODEL  PARAMETERS  FOR  SAMPLE  PAINTS 


Material 


Parameter 

A02017-001 

A02018-001 

A02018-002 

11 

1.65 

1.65 

1.65 

k 

0 

0 

0 

pxi 

_ * 

_ * 

0.044 

PX2 

_ * 

_ * 

0.044 

0 

.0064 

0.007 

0.05** 

V 

T 

15 

15 

15 

a 

40 

40 

40 

m 

1 

1 

1 

f«) 

1 

1 

1 

"a 

— 

— 

— 

\  n  n  n  n/  n 

X 

1,06  pm 

1.06  pm 

1.06  pm 

♦This  material  is  run  with  the  non- Lambertian  volume  model;  therefore 

p  values  are  not  necessary. 

X 


♦♦Material  2018-002  was  run  with  the  Lambertian  volume  model;  there 
fore  /:>v  should  not  be  used. 


TABLE  II.  TRUE  SOURCE  POLARIZATION  ANGLES 


Receiver 

Nominal  Angles 

Azimuth 

Plane 

1(0°) 

11(90°) 

+45° 

-45° 

0°-180°l 

30°-210° 

► 

5° 

-85° 

-40° 

60°-249°^ 

90°-270° 

5° 

+95° 

+50° 

Fixed  Bistatic 

5° 

+95° 

+50° 

6.1.  RE1  LECTANCE  FOR  SAMPLE  MATERIAL  A02018-001 

Material  A02018-001  is  a  green  painted  surface.  The  zero  bistatic  measurement  with  5° 
polarization  angle  (i.e.,  almost  perpendicular  polarizatton)  is  shown  in  Fig.  4.  (The  zero  bi¬ 
static  data  with  parallel-polarized  source,  although  not  shown,  have  identical  characteristics.) 
The  zero  bistatic  plot  is  sharply  peaked  at  0°,  falling  off  rapidly  to  a  constant  value  at  about 
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BIDIRECTIONAL  REFLECTANCE  (STER.’1) 


A02018  001 


FIGURE  4.  FIXED  BISTATIC  p'  FOR  A02018-001 


cos 


ti  used  in  the  model 
n 


20  .  In  all  receiver  polarizations,  p'  shows  an  angular  dependence  clearly  departing  from 

Lambertian  behavior.  The  table  of  values  for  p'  ^ 

'  n  n  n  n/ 

was  obtained  from  this  measurement  by  reading  off  pj  .and  p’  ..  at  each  angle  and  then  calcu- 

O  11  1  I* 

lating  (p’  -  p'  )  cos^  o  where  9  is  the  angle  that  the  normal  to  the  reflecting  facet  makes 

ii  ill  ft  ft 

with  the  fixed  z-axis.  In  zero  bistatic  scans,  9.  =  =  6 ^  (see  Fig.  1).  (Physically  the  source 

and  receiver  were  separated  by  1.8°.  Thus,  both  were  0.9°  from  the  true  9  .  In  the  calcula- 

n 

tions,  the  axis  was  translated  to  bring  the  x-axis  into  correspondence  with  6  =  0.)  The  sub- 

n 

traction,  pj^  -  p|  (|  eliminates  the  diffuse  contribution  which  would  distort  the  value  for 
P'(o  ,<t>  ;9  ,0  )  which  is  what  must  be  measured  (recall  Eq.  9). 


ft 

In 


ft  i 


Figs.  5,  7,  9  and  11,*  plots  of  measured  data  are  shown  for  d.  =  40,  <£.  =  180°  and  where 
0  is  scanned  in  azimuth  planes  represented  by  <f>r  =  0°,  180°;  90°,  270°;  30°,  210°;  60°,  240°. 
Each  measurement  plot  is  followed  by  plots  of  data  generated,  respectively,  by  the  Lamber¬ 
tian  model  with  no  shadowing  and  obscuration  factor,  by  the  non -Lambertian  model  with  no 
shadowing  and  obscuration  factor,  and  by  the  non- Lambertian  model  with  the  shadowing  and 
obscuration  factor.  For  example.  Fig.  6  shows  the  calculated  p'  data  for  0.  =  40°  and  6  as 
scanned  in  the  0°  and  180°  azimuth  planes  for  the  above  variations  of  the  model.  The  simu¬ 
lated  source  is  taken  to  have  a  "perpendicular"  polarization  angle.  In  these  in-plane  scans 
(4>r  =  0°,  180°),  the  main  peak  is  in  the  0°  azimuth  plane  which  is  the  forward  direction  for  the 
source  angle  of  t/>.  --  +180°.  Note  the  rise  (in  the  plot  of  measured  data)  at  large  zenith  angles 
for  the  cross-polarized  component.  This  is  a  characteristic  which  suggests  the  need  for  the 
non -Lambertian  volume  model. 


Surface  Plus  Lambertian  Volume  Model  with  No  Shadowing  and  Obscuration  Correction. 
Figure  6  plots  (in  solid  lines)  the  model  calculation  using  the  surface  model  plus  the  Lamber¬ 
tian  volume  model  with  no  correction  for  shadowing  and  obscuration.  The  following  charac¬ 
teristics  should  be  noted: 


(1)  In  the  <?>r  -  0  (forward  scattering)  azimuth  plane,  the  model  fits  the  measured  data 
very  well  between  t)  -  0  and  o  a  50°  for  matched  polarization  of  source  and  receiver. 
At  o  -  60°,  the  calculated  curve  suddenly  diverges.  This  is  thought  to  be  the  result 
of  the  failure  to  account  for  shadowing  and  obscuration  as  discussed  earlier.  At  9 

0°  and  on  into  the  backscattered  (<t>r  =  180°)  direction,  the  calculated  values  lie 
above  the  measured  values  and  this,  too,  is  believed  to  be  the  result  of  the  lack  of  a 
shadowing  and  obscuration  correction. 

(2)  In  the  cross -polarization  component  (lit),  the  model  predicts  a  flat  response  except 
for  a  slight  hump  under  the  specular  peak.  The  measured  data,  however,  show  a 
clear  angular  dependence  on  0^ 

‘Note:  On  all  reprints  of  original  computer  plots,  the  symbols  0  and  <4  are  represented 
by  and  respectively.  r 
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FIGURE  6.  CALCULATED  p'  FOR  A02018-001  USING  LAMBERTIAN  VOLUME  MODEL  AND 
NON -LAMBERTIAN  MODEL  WITH  AND  WITHOUT  SHADOWING  AND  OBSCURATION  FACTOR. 

0  =  40°,  «r  =  0°.  180°. 


22 


With  the  exception  of  these  two  characteristics,  however,  the  surface  plus  Lambertian  volume 
model  with  no  shadowing  and  obscuration  correction  fits  the  measured  data  fairly  well. 

Non- Lambertian  Volume  Model  with  No  Shadowing  and  Obscuration  Correction.  The  dotted 
lines  in  Fig.  6  show'  a  model  plot  using  the  same  parameters,  except  that  the  non- Lambertian 
volume  model  is  now  used.  Keep  in  mind  that  one  may  use  the  non -Lambertian  volume  scat¬ 
tering  as  a  model  by  itself  or  in  conjunction  with  a  specular  component.  The  latter  is  used 
here.  In  the  like-polarized  component,  nothing  has  changed  from  the  previous  case.  However, 
the  cross-polarized  component  now  fits  the  measured  data  much  more  closely.  It  rises  steadily 
at  large  angles,  both  in  the  back-scattered  and  forward-scattered  directions — a  result  of 
l/(cos  +  cos  8  )  dependence  shown  in  Eq.  (28)  for  the  volume  model.  However,  the  response 
for  the  like -polarized  component  does  not  drop  sharply  enough  at  either  side  of  the  peak,  and  at 
high  angles  in  the  forward-scattered  direction,  the  awkward  divergence  still  appears  at  60°. 

Thus,  the  non-Lambertian  volume  model  improves  the  cross -polarized  fit  (with  respect  to 
material  A02018-001)  over  that  of  the  Lambertian  volume  model  and,  apart  from  anomalies  at 
high  angles  and  near  0°,  provides  a  reasonable  fit  to  the  measurements. 

Non-Lambertian  Volume  Model  with  Shadowing  and  Obscuration  Correction.  The  dashed - 
line  curves  in  Fig.  6  show  results  with  the  shadowing  and  obscuration  correction  applied  to  the 
non-Lambertian  volume  model  calculation.  The  cross-polarized  component  is  unaffected.  The 
net  effect  on  the  match-polarized  component  is  to  reduce  the  reflectance  everywhere  except  at 
the  specular  peak  and  at  the  direct  backscattering  peak  (i.e.,  at  (i  =  0).  In  particular,  it  lowers 
the  forward-scatter  contributions  beyond  50°,  bringing  the  model  closer  to  measured  data  in 
this  region.  Overall,  the  fit  obtained  using  the  volume  model  with  a  shadowing  and  obscuration 
correction  agrees  closely  with  measurements. 

The  foregoing  discussion  applies  to  "in  plane"  receiver  scans — those  in  the  <#>r  =  0  and 
-  180°  azimuth  planes.  The  azimuth  plane  perpendicular  to  the  0°,  180°  plane  is  the  90°, 

270°  plane  and  is  referred  to  as  "out-of-plane".  The  plane  we  are  in  or  out  of  is  the  plane  of 
incident  beam  and  target  normal,  or  the  target  incidence  plane.  (See  Fig.  1.) 

In  Fig.  7  we  have  the  plot  of  measured  data  for  the  out -of -plane  situation  with  perpendicular  - 
polarized  source  again.  In  this  case,  however,  the  incidence  plane  is  perpendicular  to  the  re¬ 
flection  plane.  At  o,  therefore,  plane  is  the  same  as  out  of  plane.  For  this 

reason,  the  reflectances  of  match-polarized  and  cross -polarized  components  seem  to  exchange 
behaviors  in  the  out-of-plane  configi ration,  as  is  verified  by  the  plotted  measurements  as  well 
as  by  the  model  calculations.  Figure  8  presents  plots  of  a  Lambertian  model  without  the  shad¬ 
owing  and  obscuration  factor,  a  surface  plus  non-Lambertian  volume  model  without  the  shad¬ 
owing  and  obscuration  factor,  and  the  surface  plus  non-Lambertian  volume  model  with  the 
shadowing  and  obscuration  factor.  As  before,  it  is  apparent  that  the  use  of  the  non-Lambertian 
volume  model  plus  the  shadowing  and  obscuration  factor  improves  agreement  between  model 
and  measurements  so  that,  apart  from  a  possible  overall  scale  factor,  the  agreement  is  within 
measurement  fluctuation.  23 
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FIGURE  7.  MEASURED  p’  FOR  A02018-001.  =  40°,  «  =  180°,  <t>  r  =  90°,  270°. 
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FIGURE  8.  CALCULATED  p'  FOR  A02018-001  USING  LAMBERTIAN  VOLUME  MODEL  AND 
NON- LAMBERTIAN  VOLUME  MODEL  WITH  AND  WITHOUT  SHADOWING  AND  OBSCURA¬ 
TION  FACTOR.  -  40°,  =  180°;  *  =  90°,  270°. 
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For  additional  validation,  plots  are  shown  for  the  30°,  210°  azimuth  planes  (Figs.  9  and  10) 
and  for  the  60°,  240°  azimuth  planes  (Figs.  11  and  12).  The  characteristics  of  calculated  and 
measured  curves,  apart  from  a  scale  factor,  are  in  excellent  agreement.  Figures  13  through 
20  represent  similar  comparisons  for  the  case  when  the  source  polarizer  is  set  for  -45°  (in 
the  0°,  180°  azimuth  plane)  and  set  parallel  (in  the  30°,  210°;  60°,  240°;  and  90°,  270°  planes). 
Measured  plots  are  presented  with  the  calculated  plot  to  represent  the  surface  plus  non- 
Lambertian  volume  model  and  to  include  the  shadowing  and  obscuration  factor. 

6.2.  REFLECTANCE  FOR  SAMPLE  MATERIAL  A02018-002 

Material  A02018-002  consists  of  a  tan  painted  surface. 

Based  on  the  zero  bistatic  scan,  Fig.  21,  material  A02018-002  appears  to  be  somewhat 
brighter  than  material  A02018-001.  Whereas  the  non-Lambertian  volume  model  was  clearly 
the  best  choice  for  material  A020 18-001,  it  is  not  in  the  case  of  A02018-002.  In  this  latter 
case,  the  best  choice  is  the  Lambertian  model. 

The  lack  of  angular  dependence  in  the  reflectance  of  the  cross -polarized  component  could 
have  a  number  of  explanations.  Multiple  scattering  increases  for  rougher  surfaces.  Since  such 
scattering  may  not  be  angular  dependent,  it  could  become  a  large  enough  factor  to  swamp  the 
angular  dependence  which  is  otherwise  present.  Moreover,  the  difference  in  color  between  the 
green  and  tan  certainly  alters  the  absorption  and,  consequently,  can  alter  the  angular  dependence 
as  well. 

In  any  case,  the  appropriate  model  to  use  can  be  determined  by  looking  at  the  cross- 
polarized  component  of  the  fixed  bistatic  scan.  If  a  clear  angular  dependence  is  present,  the 
non-Lambertian  model  should  be  used.  But  if  there  is  little  or  no  apparent  angular  dependence, 
as  with  material  A02018-002,  then  the  Lambertian  model  is  more  appropriate. 

In  Figs.  22  through  29,  plots  are  provided  for  different  azimuth  planes,  beginning  with  the 
plot  for  measured  data,  followed  immediately  by  the  corresponding  plot  from  model  calcula¬ 
tions.  In  this  group  of  illustrations,  Figs.  22  through  25  represent  perpendicular  source  po¬ 
larization,  while  Figs.  26  through  20  represent  a  source  parallel  polarization  for  the  0°,  180° 
azimuth  plane  and  for  the  90°,  270°  azimuth  plane. 

In  all  cases  the  fit  appears  to  be  excellent,  except  for  occasional  anomalies  at  large  azi¬ 
muth  angles.  Further  modification  of  the  shadowing  and  obscuration  factor  should  decrease 
these  present  anomalies. 

6.3.  POLARIZATION  ANGLE  (^f)FOR  SAMPLE  MATERIAL  A020I8-001 

The  reflectances  of  the  perpendicular  and  parallel  components  of  a  linearly  polarized 
beam  vary  as  functions  of  the  source-receiver  angles  and  the  index  of  refraction  of  the  target 
material.  (See,  for  example,  the  Fresnel  equations,  Ref.  3.)  Based  on  observations,  the  index 
of  refraction  varies  little  over  a  wide  range  of  paint  surfaces.  For  the  particular  materials 
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FIGURE  9.  MEASURED  p'  FOR  A02018-001.  6  =  40°, 


180°,  *  -  30°,  210°. 
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FIGURE  10.  CALCULATED  p'  FOR  A02018-001  USING  NON -LAMBERTIAN  vOLUME  MODEL 
WITH  SHADOWING  AND  OBSCURATION  FACTOR.  =  40°.  ^  =  180°,  «r  -  30°,  210°. 
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FIGURE  17.  MEASURED  p'  FOR  A02018-001.  9{  =  40°,  ^  =  180°,  <t>f  *  30,  210°. 
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FIGURE  18.  CALCULATED  p'  FOR  A02018-00I  USING  NON -LAMBERTIAN  VOLUME  MODEL 
WITH  SHADOWING  AND  OBSCURATION  FACTOR.  0[  =  40°,  *  =  180°,  *  =  30°,  210°. 
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FIGURE  21.  FIXED-BISTATIC  p’  FOR  A02018-002 
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FIGURE  22.  MEASURED  p'  FOR  A02018-002.  ^  =  40°,  ^  =  180°,  =  0°,  180°. 
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FIGURE  23.  CALCULATED  p’  FOR  A02018-002  USING  LAMBERTIAN  VOLUME  MODEL. 

S,  »  40°,  «  180°,  =  0°,  180°. 
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FIGURE  25.  CALCULATED  p'  FOR  A02018-002  USING  LAMBERTIAN  VOLUME  MODEL. 

=  40°,  0j  =  180°,  0r  =  90°,  270°. 
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FIGURE  26.  MEASURED  p'  FOR  A02018-002.  0.  --  40°,  *  =  180°.  <*>r  =  0°,  180°. 
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FIGURE  28.  MEASURED  p'  FOR  A02018-002.  d{  =  40°,  =  180°,  4>r  =  90°,  270°. 
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covered  in  this  report,  the  variation  is  considered  to  be  zero.  Therefore,  for  these  surfaces, 
polarization  angle  is  essentially  a  function  only  of  source-receiver  positions. 

As  additional  validation  for  the  model,  predicted  polarization  angles  are  compared  with 
polarization  angles  extracted  from  the  measured  data.  Figures  30  through  33  show  plots 
obtained  for  the  0°,  180°;  90°,  270°;  30°,  210°;  and  60°,  240°  azimuth  planes.  Measured  data 
represent  material  A02018-001.  In  all  cases,  agreement  between  measurements  and  calcu¬ 
lations  is  excellent,  with  the  average  disparity  not  more  than  10%.  In  particular,  the  dramatic 
agreement  between  measurements  and  model  in  the  30°,  210°  and  60°,  240°  azimuth  planes 
constitutes  powerful  verification  of  the  model  and  affirms  its  usefulness  in  arbitrary  source- 
receiver  positions. 

6.4.  PERCENT  POLARIZATION  FOR  SAMPLE  MATERIALS  A02018-001  AND  A02018-002 

Percent  polarization  (P)  validates  the  ratio  of  surface -to -volume  contributions  to  reflec¬ 
tance.  Percent  polarization  depends  on  both  polarized  reflectance  and  angle  of  polarization, 
both  validated  in  earlier  sections  of  this  report.  In  this  section,  we  compare  model  predic¬ 
tions  with  percent  polarization  values  extracted  from  measured  data. 

Figures  34  and  35  illustrate  degree  of  polarization  for  scans  of  material  A02018-002,  for 
perpendicular  and  parallel  sources,  respectively.  The  validity  of  the  model  is  supported  by 
the  close  correlation  between  the  behavior  of  values  extracted  from  measured  data  and  those 
calculated  with  the  model. 

Additional  confirmation  of  the  model  is  provided  in  Figs.  36  through  38  where  percent 
polarization  plots  are  given  for  material  A02018-001  in  the  0°,  180°  and  90°,  270°  azimuth 
angle  planes. 

6.5.  REFLECTANCE  FOR  SAMPLE  MATERIAL  A02100 

Here  the  material  was  a  soil  specimen.  The  fixed -bistatic  scan  (see  Fig.  39)  indicates  a 

strong  angular  dependence  in  both  the  like -polarized  and  the  cross-polarized  components  and 

no  specularity.  Moreover,  the  angular  dependence  in  the  fixed -bistatic  scan  looks  very  much 

like  the  — - - -r-  dependence  in  the  non -Lambertian  volume  model.  However,  in  the  0 

cos  9.  +  cos  9  r 

I  r 

scans  with  9^  held  fixed,  the  angular  dependence  is  no  longer  typical  (see  Figs.  40,  43,  and  46). 

Since  there  is  no  apparent  specularity  in  the  measurement  data,  the  model  was  run  so  as 
to  consider  only  Lambertian  or  non -Lambertian  components,  with  no  specular  component.  Our 
validation  was  done  with  a  perpendicular -polarized  source  for  =  0°,  20°,  and  40°,  respectively, 
and  #r  scanned  in-plane.  For  each  0y  the  measurement  graph  is  given  first,  followed  by  the 
graphs  of  the  non -Lambertian  model  and  the  Lambertian  model.  Note  that  in  selecting  p 
value  for  the  Lambertian  model,  we  must  take  an  average  of  the  cross -polarized  part  of  the 
fixed -bistatic.  Therefore,  this  value  is  slightly  higher  than  the  py  used  in  the  non -Lambertian 
model.  (Section  7  describes  how  to  select  parameters.) 
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For  the  0.  =  0°  case,  note  that  the  measured  data  (see  Fig.  40)  lacks  the  characteristic 
angular  dependence  of  the  non -Lambertian  model  (see  Fig.  41)  —  in  fact,  the  graph  falls  off 
below  the  Lambertian  graph  (see  Fig.  42). 

For  the  &  =  20°  case,  the  measured  data  (see  Fig.  43)  shows  an  upturn  at  higher  0 values, 
which  is  more  consistent  with  the  non-Lambertian  model.  In  fact,  the  non -Lambertian  model 
prediction  (see  Fig.  44)  shows  much  closer  agreement  than  in  the  0.  =  0°  case.  The  Lambertian 
calculation  (see  Fig.  45)  once  again  appears  to  be  taking  a  rough  average.  The  situation  for 
0.  =  40°  (see  Figs.  46-48)  is  very  similar  to  that  for  0.  =  20°. 

Behavior  of  the  fixed -bistatic  measurement  data  indicates  that  the  non-Lambertian  compo¬ 
nent  should  dominate.  However,  at  small  values  of  0.  the  accuracy  is  not  good.  It  is  not  abso¬ 
lutely  clear  where  the  difficulty  lies.  One  should  note,  however,  that  we  took  the  depolarization 
to  be  identically  equal  to  1.  Further  modeling  to  determine  the  depolarization  dependence  more 
accurately  may  well  resolve  the  problem. 
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FIGURE  30.  VARIATION  OF  POLARIZATION  ANGLE  OF  REFLECTED  RADIANCE  AS 
FUNCTION  OF  SOURCE -RECEIVER  POSITION.  8  =  40°,  *  =  180°,  <t>r  =  0°,  1800. 
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▼ :  Measured  Source,  -40 
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FIGURE  31.  VARIATION  OF  POLARIZATION  ANGLE  OF  REFLECTED  RADIANCE  AS 
FUNCTION  OF  SOURCE -RECEIVER  POSITION.  =  40°,  =  180°,  *r  =  30°.  210®. 


O :  Measured  Source,  +5° 
V:  Measured  Source,  -40° 
O :  Measured  Source,  -85° 
x :  Calculated  Source,  +5° 
©:  Calculated  Source,  -40° 
▲  :  Calculated  Source,  -85° 
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FIGURE  33.  VARIATION  OF  POLARIZATION  ANGLE  OF  REFLECTED  RADIANCE  AS 
FUNCTION  OF  SOURCE -RECEIVER  P06ITI0N.  »  40°,  *  s  180°,  *r  =  90°,  270°. 


FIGURE  34.  PERCENT  POLARIZATION 
VARIATION  FOR  A02018-002  AS  FUNC¬ 
TION  OF  SOURCE -RECEIVER  POSITION. 
Q  =  40°,  *  =  180°,  <t>r  =  0°,  180°;  perpen- 
1  dlcular  source. 


FIGURE  35.  PERCENT  POLARIZATION 
VARIATION  FOR  A02018-002  AS  FUNC¬ 
TION  OF  SOURCE -RECEIVER  POSITION. 
=  40°,  =  180°,  <t>r  =  0°,  180°;  parallel 

source. 


VARIATION  FOR  A020I8-001  AS  FUNC¬ 
TION  OF  SOURCE -RECEIVER  POSITION. 
3  40°,  =  180°,  4r  =>  0°,  180°;  perpen¬ 

dicular  source. 
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FIGURE  37.  PERCENT  POLARIZATION 
VARIATION  FOR  A02018-001  AS  FUNC¬ 
TION  OF  SOURCE -RECEIVER  POSITION. 
e  =  400,  0  =  180°,  0  =  90°,  270°. 


FIGURE  38.  PERCENT  POLARIZATION 
VARIATION  FOR  A020I8-001  AS  FUNC¬ 
TION  OF  SOURCE -RECEIVER  POSITION. 

0.  =  40°,  0.  -  1800,  0  90°,  270°. 
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BIDIRECTIONAL  REFLECTANCE  CSTER.‘») 


FIGURE  39.  FIXED  BISTATIC  p'  FOR  A02100 


BIDIRECTIONAL  REFLECTANCE  (STER."!) 


R02100  101 


FIGURE  40.  MEASURED  p'  FOR  A02I00.  9{  =  0°,  ^  =  0°,  =  0°,  180°. 
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BIDIRECTIONAL  REFLECTANCE  (STER.-‘I 


R02100  101 


FIGIRE43.  MEASURED FOR  AU2  .*•>  20°.  .■  1*>0U  U|,  0°,  I8I)J 


BIDIRECTIONAL  REFLECTANCE  (STER.“l) 


FIGURE  44.  CALCULATED  p'  FOR  A02100  USING  NON -LAMBERTIAN  VOLUME 
MODEL,  d.  =  20°,  *  =  180°,  <t>f  =  0°,  180°. 
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BIDIRECTIONAL  REFLECTANCE  (STER.*1) 
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FIGURE  46.  MEASURED  p1  FOR  AOiGOO,  ti  -  40°,  <*.  =  180°,  *r 
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BIOIRECTIONHL  REFLECT fiNCE  (STEP.”1) 
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FIGURE  47.  CALCULATED/)'  FOP  A02100  USING  NON -LAMBERTIAN  VOLUME 
MODEL.  •*.  40°,  <i>.  180°,  </>  ,  0°,  180°. 


BIDIRECTIONAL  REFLECTANCE  (STER.'1) 


FIGURE  48.  CALCULATED  p'  FOR  A02100  USING  LAMBERTIAN  VOLUME  MODEL. 

8,  -  40°.  <b,  -  180°,  «  =  0°.  180°. 
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MODEL  PARAMETERS 


This  section  briefly  describes  the  model  ,>.i.  imeters  that  can  be  used  in  the  bidirectional 
reflectance  program  and  explains  how  their  values  are  derived.  The  choice  of  parameters  for 
use  in  the  program  depends  to  some  extent  on  the  mode  of  the  model  being  run.  Basically,  the 
model  is  run  in  three  different  modes: 

(1)  Surface  and  Lambertian  volume  components 

(2)  Aou -Lambertian  volume  component  too  surface  contribution  included) 

(3)  Surface  and  non-Lambert  inn  volume  ,  ompouents 

Therefore,  we  have  grouped  the  model  parameters  as  follows: 

(1)  Polarization  parameters 

(2)  Surface  model  parameters 

(3)  Lambertian  volume  model  parameter. 

(4)  Non-Lambertian  volume  model  parameters 

7.1.  SOURCE  POLARIZATION  PARAMETERS 

The  present  model  has  been  designed  to  account  for  polarization  dependence  in  both  sur¬ 
face  and  volume  components. 

In  the  suriace  component,  polarization  is  accounted  lot  automatically  in  the  Fresnel  reflec¬ 
tance  coefficients.  In  the  most  general  case,  .  uen  polarization  can  be  elliptical  and  can  be  de¬ 
composed  into  linear  and  circular  components.  To  date,  only  a  linearly  polarized  source  and 
receiver  have  been  used  in  measurements.  However,  for  some  applications,  circularly  polar¬ 
ized  sources  or  receivers  may  be  of  interest.  Therefore,  in  the  model,  we  have  provided  pro¬ 
gram  subroutines  which  take  into  consideration  the  eltipticity  and  handedness  (i.e,,  direction 
of  rotation  in  an  elliptically  polarized  source)  of  noth  incident  ana  retlected  beam. 

For  volume  components  in  both  Lamoertiaa  and  iion-Lambeman  cases,  it  is  assumed  that 
'ellectam  e  will  be  depolarized  to  some  extern.  In  butn  cases,  in  fact,  we  assume  total  depolar¬ 
ization.  Therefore,  although  a  depolarization  iacior  has  been  included  in  the  non -Lambertian 
volume  model  for  future  flexibility,  we  assume  DP(,D  1. 

The  source  polarization  may  most  generaii;  be  defined  as  partially  polarized  with  the  po¬ 
larized  component  elliptically  polarized.  The  state  of  polarization  of  tne  source  will  be  defined 
by  its  degree  of  polarization,  P.  and  parameters  A,  B.  ^ .  and  H  io  detine  the  elliptical  polariza¬ 
tion  of  the  polarized  component.  Here.  A  and  B  are  the  intensities  along  the  semi -major  and 
semi-minor  axes,  respectively.  The  angle  v  ^  the  angle  the  semi-major  axis  of  the 

ellipse  and  tne  direction  normal  to  the  plane  m  incidence,  measured  looking  into  the  source 
beam;  ^  is  equivalent  too  except  that  0°  •  v  »th>°  and  -tfOu  ..  9uu  1'iie  handedness  H 
-  *1. 


do 


The  Stokes  vectors  provide  a  convenient  formalism  for  defining  the  polarization  state  of 
the  reflected  radiance.  (Reference  [5]  provides  a  general  discussion  of  Stokes  vectors  in  this 
context.  ) 


Ip  cos  2x  cos  2i p 

I  cos  2x  sin  2<^ 
P 

I  sin  2x 
P 


where  I  and  I  are  the  polarized  and  unpolarized  components,  respectively,  in  the  reflected 
p  u 

radiance.  The  degree  of  polarization  in  the  reflected  radiance  is  P  =  1/(1  +  I  ).  Angles  x 

r  r  u 

and  define  the  polarization  state  of  the  reflected  radiance:  is  the  angle  between  the  semi¬ 

major  axis  of  the  ellipse  and  the  direction  normal  to  the  plane  of  reflection;  tan  x  =  ±Vb/ A 
where  A  and  B  are  the  intensities  along  the  semi-major  and  semi -minor  axes  of  the  polariza¬ 
tion  ellipse  and  tan  x  >  q  for  |®^"t_handed  elliptically  polarized  radiation. 

The  RHOPRIME  program  produces  the  Stokes  vector  S  for  unit  irradiance  in  the  input  beam; 
the  area  may  also  be  defined  to  be  unity  and  then  S  represents  a  reflectance  Stokes  vector.  The 
program  also  produces  the  components  of  the  reflected  radiance  transmitted  with  a  receiver 
polarization  analyzer  oriented  parallel  or  perpendicular  to  the  reflectance  plane  for  computing 

'Wdp*,r 


7.2.  SURFACE  MODEL  PARAMETERS 

n  and  k.  These  are  the  real  and  imaginary  parts  of  the  refractive  index.  As  discussed 
earlier  in  this  report,  they  are  used  for  the  determination  of  R(/i),  the  Fresnel  reflectance. 
Values  for  n  and  k  are  estimated  for  the  paint  surfaces  *n  this  study.  Moreover,  the  surface 
is  assumed  to  be  essentially  nonconducting  so  that  k  =  0-  Based  on  experience  with  similar 
paint  samples,  n  is  taken  to  be  1.65.  For  a  given  sample,  n  and  k  can  be  determined  accu¬ 
rately  by  measuring  the  Brewster  angle  and  calculating  n  and  k  as  outlined  in  Section  4  on  the 
surface  model.  At  the  present  time  the  program  used,  RHOPRIME,  does  not  do  this. 

t  and  fl .  These  parameters  are  used  in  the  function  which  provides  a  correction  to  the 
program  to  account  for  shadowing  and  obscuration  effects  result lr  ^  from  the  roughness  of  the 
surface.  Values  for  r  and  Q  have  been  selected,  base  1  on  observed  characteristics  of  reflec¬ 
tance  properties.  They  have  been  established  as  r  -  15  and  =  40. 

p'  fa  ,$  ;  ,$  \c°s2  o  •  One  of  the  quantities  in  Eq.  (9)  from  which  e  ,4  )  is  de- 

\ft  n  ft  ft/  ft  ‘  ’  r  r 

termlned  is  p'/e  ,$  :  \cos^  As  previously  discussed,  p'/e  ,<*».;  (<A,$  \  is  obtained  from 

\  ft  ft  n  h  J  n  2  \nnrin/ 

from  zero  bistatlc  data.  Values  for  p'/h  ,$  ;  v  ,<p  \cos^  o  must  be  calculated  (preferably  for 

I  A  AAA  A 

\  n  n  n  nl  n 

increments  of  two  degrees)  and  made  into  a  table  which  is  one  of  the  mod.  I  inputs.  Values  are 
provided  in  Tables  lll-VI  for  materials  A02017-001,  A02U18-001,  A02018-002,  and  A02100 
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7.3.  LAMBERTIAN  VOLUME  MODEL  PARAMETERS 


and  These  are  the  cross  components  of  polarized  radiation  used  in  the  model  to 

account  for  the  diffuse  contribution,  p^j  2p(| .  and  p^9  2 p^  ()  ,  where  pjj  and/or  p|(^  is  de¬ 

termined  by  taking  the  average  value  of  the  ci  „ss  component  from  the  measured  data.  Ac¬ 
cording  to  the  reciprocity  theorem,  pAl  p^  j  Ref.  ti| .  It  is  important  to  remember  that  p^ 
values  only  used  when  the  volume  scatter  model  is  not  used.  When  the  volume  model  is  used, 


7.4.  NON- LAMBERTIAN  VOLUME  MODEL  PARAMETERS 

p  This  represents  the  non- Lambertian  volume  scatter  component;  it  is  determined  by 
extracting  p'  ^  or  p’  at  the  point  which  would  lie  under  the  peak  of  the  zero  bistatic  scan  if  the 
measured  curve  were  smooth;  p^  -  2pjn  -  2p  ( ,  at  the  peak  point.  (The  fact  that  a  hump  some¬ 
times  occurs  on  the  measured  cross  component  curve  is  discussed  in  the  section  on  Model  Vali¬ 
dations.) 

Here  again,  it  is  important  to  remember  that  when  the  Lambertian  volume  model  is  used, 
p  -  0  and  p  *  0.  When  the  non- Lambertian  volume  model  is  used,  p  .  /  0  and  p  =-  0.  Also, 

p  and  p  are  never  simultaneously  nonzero  m  models  which  have  been  validated  to  date. 

V  X 

DP(Pi,  f(,i),  Y  Integral  parts  of  SUBROUTINE  FUNC,  these  parameters  currently  are 
ail  set  equal  to  1.  They  have  been  included  to  provide  flexibility  for  later  model  modifications. 


table  in.  p'(0n,tf>n;  0n,*n)  cos2  #n  values  for  A02017-001 


e. 

_ n 

P’  ^6n»^n»6«»^^COS 

n  n  n  n 

.9 

.14528 

2.9 

.08415 

4.9 

.0529 

6.9 

.03705 

8.9 

.03077 

10.9 

.02315 

12.9 

.01829 

14.9 

.01786 

16.9 

.01524 

18.9 

.01344 

20.9 

.01169 

22.9 

.01091 

24.9 

.01148 

26.9 

.00984 

28.9 

.0086 

30.9 

.00925 

32.9 

.00908 

34.9 

.00815 

36.9 

.00625 

38.9 

.00649 

40.9 

.00664 

42.9 

.00617 

44.9 

.00531 

46.9 

.00465 

48.9 

.00474 

50.9 

.00406 

52.9 

.00354 

54.9 

.00291 

56.9 

.00323 

58.9 

.00265 

60,9 

.00224 

62.9 

.00217 

64.9 

.00211 

66.9 

.00205 

66.9 

.00196 

70.9 

.0027 

69 


TABLE  IV.  p'(0n,4>n;  0n,$n)  cos2  0n  VALUES  FOR  A02018-001 


e- 
_ n 

p* 

n  n  n  n 

1.1 

.5143 

3.1 

.27294 

5.1 

.11897 

7.1 

.05797 

9.1 

.03291 

11.1 

.02348 

13.1 

.02058 

15.1 

.01614 

17.1 

.0143 

19.1 

.0118 

21.1 

.01191 

23.1 

.01056 

25.1 

.01055 

27.1 

.00924 

29.1 

.00918 

31.1 

.00672 

33.1 

.00708 

35.1 

.00612 

37.1 

.00597 

39.1 

.00578 

41.1 

.00549 

43.1 

.00469 

45.1 

.00459 

47.1 

.00433 

49.1 

.00417 

51.1 

.0046 

53.1 

.00376 

55.1 

.00292 

57.1 

.0023 

59.1 

.00179 

61.1 

.00168 

63.1 

.00132 

65.1 

.0013 

67.1 

.00109 
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TABLE  V.  p’(0  ,<f>  ;  6  ,<fr  )  cos2  6  VALUES  FOR  A02018-002 
H  n’*n  n’*n  n 

°A  p’(6-,4.-;e-,»-)Cos28- 


1.15 

.02497 

3.15 

.02520 

5.15 

.02388 

7.15 

.01975 

9.15 

.02135 

11.15 

.02177 

13.15 

.02335 

15.15 

.02353 

17.15 

.02300 

19.15 

.02243 

21.15 

.02338 

23.15 

.02022 

25.15 

.02056 

27.15 

.01726 

29.15 

.01885 

31.15 

.01966 

33.15 

.01875 

35.15 

.02089 

37.15 

.01882 

39.15 

.01709 

41.15 

.01785 

43.15 

.01604 

45.15 

.01545 

47,15 

.01295 

49.15 

.01319 

51.15 

.01434 

53.15 

.01180 

55.15 

.01385 

57.15 

.01190 

59.15 

.00922 

61.15 

.01149 

63.15 

.00869 

65.15 

.00954 

67.15 

.00808 

69.15 

.01349 
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TABLE  VI.  P'(0n,0n;  0n,*n)  cos2  ^  VALUES  FOR  A02100 


e 

_ n 

p,(0n»*n;On’*n 

0.1 

.21399 

2.1 

.16918 

4.1 

.1096 

6.1 

.08454 

8.1 

.05494 

10.1 

.03917 

12.1 

.03869 

14.1 

.02547 

16.1 

.03465 

18.1 

.02112 

20.1 

.02753 

22.1 

.02155 

24.1 

.02261 

26.1 

.01906 

28.1 

.02055 

30.1 

.01637 

32.1 

.01799 

34.1 

.01778 

36.1 

.01601 

38.1 

.01559 

40.1 

.01337 

42.1 

.01626 

44.1 

.01101 

46.1 

.01054 

48.1 

.00894 

50.1 

.01003 

52.1 

.00898 

54.1 

.00648 

56.1 

.00815 

58.1 

.00701 

60.1 

.00470 

62.1 

.00552 

64.1 

.00497 

66.1 

.00387 

68.1 

.00343 

70.1 

.00292 
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Appendix  I 

DOCUMENTATION  OF  BIDIRECTIONAL  REFLECTANCE  PROCRAM  (RHOPRIME) 

Program  RHOPRIME  is  the  main  calling  program  for  subroutines  to  read  and  store  ma¬ 
terials  data,  perform  geometrical  calculations,  compute  bidirectional  reflectances  for  any 
source/receiver  position  and  polarization,  and  prepare  the  output  in  a  convenient  format.  The 
calling  sequence,  purpose,  and  calculations  performed  by  each  subroutine  are  given  below,  fol¬ 
lowed  by  details  on  the  input  data  formats. 


III.  1.  DESCRIPTIONS  OF  SUBROUTINES 

SUBROUTINE  INDATA.  This  is  the  first  subroutine  called.  Material  parameters  needed 
for  the  calculation  of  bidirectional  reflectance  are  read  and  stored.  Material  parameters  are 


MAT  =  material  specifier 

N  =  n  =  real  part  of  refractive  index 

K  =  k  -  imaginary  part  of  refractive  index 

RX1  -  p  .  -  diffuse  reflectance  for  i  polarized  source 

A  * 

RX2  =  p  0  =  diffuse  reflectance  for  II  polarized  source 


RHOV  = 

SIGMA 

RPO 


X2 

p  =  volume  reflectance 
Kv 


l  parameters  available  to  calculate  p' It)  ,<t>  ;  0  ,<t>  \  in  subroutine  FUNC 

I  A  A  A  A  I 

J  \  n  n  n  n/ 


TAU  =  r(deg) 
OMEGA  =  Q  (deg) 
Qi 
Q2 


parameters  to  calculate  a  shadowing  and  obscuration  factor  to  be  ap¬ 
plied  to  p'  (cos  f’NP)  in  subroutine  FUNC 


RCOSBNP  -- p'/o  r  .  ;  \  COS2  V,  j 

\  ri  n  ri  n/  n  l  table  of  zero-degree  bistatic  bidirectional  ref lec- 

BNP  =  6 ^  (deg)  J  tance  data 

n  * 

TITLE.  A  title  card  (optional)  is  read  and  used  to  identify  on  the  printed  output  the  cal¬ 
culations  to  be  performed. 


FACET.  The  source  and  receiver  are  located  in  an  earth-fixed,  right-handed  XYZ  coordi¬ 
nate  system.  The  XYZ  components  of  the  unit  normal  vector  of  the  reflecting  surface  are  read 
(optional).  If  the  facet  definition  card  is  not  supplied,  the  facet  unit  normal  vector  defaults  to 
(0,  0,  1). 


COMPUTATION  REQUEST.  The  specification  of  source  and  receiver  positions  and  source 
polarization  for  computation  of  the  bidirectional  reflectance  is  read. 
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ISW  model  select  or 

TS  =  zenith  angle  of  source  (deg) 

PS  -  azimuth  angle  of  source  (deg) 

TD  -  zenith  angle  of  receiver  (deg) 

PD  -  azimuth  angle  of  receiver  (deg) 

A  -  intensity  of  major  axis  of  polarization  ellipse 
B  --  intensity  of  minor  axis  of  polarization  ellipse 

PSI  angle  of  major  axis  of  polarization  ellipse  from  the  normal  to  the  plane  of  incidence 
measured  CCW  looking  into  the  source,  0  PSI  -  180  (deg) 

P  -  polarization  of  source  (0  P  1.0) 

II  ;  handedness  of  polarization  ellipse  (±1.0  or  0.0) 

MI  material  specifier 

SUBROUTINE  SCAN.  This  subroutine  defines  a  sequence  of  detec* or  positions  for  a  speci¬ 
fied  source  position  and  polarization. 


ISW  model  selector 
TS  zenith  angle  of  source  (deg) 

PS  azimuth  angle  of  source  (deg) 

'IDS  start  zenith  angle  of  receiver  (aeg; 

TDK  end  zenith  angle  of  receiver  tdeg) 

TSTEP  zenith  angle  scan  increment  (deg) 

PDS  start  azimuth  angle  of  receiver  (cieg, 

PDL  end  azimuth  angle  of  receiver  (de 
PSTEP  azimuth  angle  scan  increment  meg; 

»A  intensity  oi  major  ;ixis  of  polarization  ellipse 
B  intensity  of  minor  axis  of  polarization  ellipse 

PSI  angle  of  major  axis  of  polarization  empse  uom  the  normal  m  me  plane  of  incidence 
measured  CCW  looking  into  me  souzet  0  PSI  130  (deg) 

P  polarization  oi  source  (0  p  1.0) 

!l  handedness  ot  polarization  ellipse  t  *  J.O  or  u 
Ml  material  speedier 

SUBROUTINE  GEOM.  Tins  subroutine  does  the  necessary  geometrical  calculations  of 
angles  needed  for  the  bidirectional  reflectance  calculations  (see  Ftg.4‘J  /. 


OR  (t),  0.  1)  is  a  unit  vector  along  the  eat  th -fixed  Z  axis 

PSI  the  angle  of  me  major  axis  of  puUn/.dion  ellipse  from  the  normal  vector  of  the  OR, 
E  plain  measured  CCW  looking  into  •  te  source,  0  PSI  180  uleg) 


X 


D  E 
D  K 
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FIGURE  49.  BIDIRECTIONAL  REFLECTANCE  GEOMETRY 
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XA 
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U  =■ 


L>  x  E 
D  x  El 


XAP 

YAP 


OP  x  E 
I  OP  *El 

OP  x  D 

[op  x  dI 


COSb  X  •  D 
COSBDP  -  OP  •  D 
COSBEP  -  OP  E 
COSBNP  -  OP  •  X 

PSIPE  -  PSI  -  SIGN(-XAP  •  OR)ARCOS(XAP  •  Y) 

-  angle  of  major  axis  of  polarization  ellipse  from  the  normal  vector  of  the  OP,  E 
plane 

PSIDE  -  PSI  -  SIGN(Y  •  D)ARCOS(U  •  Y) 

angle  of  major  axis  of  polarization  ellipse  from  the  normal  vector  of  the  D,  E 
plane 

WADE  -  -SIGN(-YA  •  E)ARCOS(XA  •  YA) 

-  angle  for  transforming  the  output  angle  of  polarization  from  E,  D  plane  to  OR,  D 
plane 

EDPH1  -  ARCOS(XAP  •  YAP)  the  relative  azimuth  angle  between  E  and  D  in  the  facet 
coordinate  system 

DC  (-S1NBEP,  0,  COSBEP)  direction  of  specular  ray  in  the  facet  coordinate  system 
Dl  (SINBDP  COS  EDPHI,  SINBDP  SIN  EDPHI,  C06BDP)  =■  direction  of  reflected  ray  in 
the  facet  coordinate  system 
NZi  DC  OP 
NZ  NZI  DC 
DN  Dl • NZ 

PHIEN  0  IF  DN  0  |_  parameter  required  in  FUNCTION  FUNC 

..  2  -  ARCOS(-DN)  IE  DN  -  fur  shadowing  and  obscuration 

SUBROUTINE  GFRM.  GFRM  does  all  of  me  bidirectional  reflectance  calculations.  1116 
surbroutinc  requires: 


F  series  of  switches  whicli  can  be  set  (optional)  to  reduce  the  number  of  redundant  compu¬ 
tations  when  GFRM  is  used  as  part  of  the  multifaceted  target  model 
COSH  defined  in  SUBROUTINE  GEOM 
COSRDP  defined  m  SUBROUTINE  GEOM 
COSBEP  defined  in  SUBROUTINE  GEOM 
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‘*S*,*,\  v 


COSBNP  -  defmed  in  SUBROUTINE  GEOM 
PSIPE  -  defined  in  SUBROUTINE  GEOM 
PSIDE  =  defined  in  SUBROUTINE  GEOM 
WADE  *  defined  in  SUBROUTINE  GEOM 

AP  =  area  of  facet  (if  AP  -  zero,  GFRM  returns  a  bidirectional  reflectance  Stokes  vector; 
if  AP  *  0,  GFRM  i  eturns  a  Stokes  vector  for  the  reflected  radiant  intensity  for  unity 
irradiance  in  the  incident  beam) 

MI  =  material  specifier  (available  in  COMMON) 

ISW  =  model  selector  (available  in  COMMON) 

W  =  wavelength  specifier  (pvailable  in  COMMON),  not  used 

TABLE  =  array  containing  idl  of  the  materials  properties  data  read  in  SUBROUTINE  IN¬ 
DATA 

GFRM  returns  the  bidirectional  reflectance  Stokes  vector  (AP  =  0)  or  radiant  intensity  Stokes 
vector  (AP  *  0). 

Ill  =  Stokes  vector  for  surface  plus  Lambertian  model  with  polarized  source 
12 1  =  Stokes  vector  for  surface  plus  Lambertian  model  with  unpolarized  source 

113  =  Stokes  vector  for  non -Lambertian  volume  model  with  polarized  source 

123  =  Stokes  vector  for  non-Lambertian  volume  model  with  unpolarized  source 

114  =  Stokes  vector  for  combined  model  with  polarized  source 

124  -  Stokes  vector  for  combined  model  with  unpolarized  source 

FUNCTION  GETDAT  returns  the  appropriate  material  parameters  for  bidirectional  re¬ 
flectance  calculations,  namely  N,  K,  RX1,  RX2,  RHOV,  RCOSBNP,  DPO,  DP90,  F,  G. 

FUNCTION  FUNC  provides  the  optional  capability  for  deriving  RCOSBNP  analytically  (if 
SIGMA  *  0)and  for  deriving  a  shadowing  and  obscuration  correction  factor  (optional)  to  the 
RCOSBNP  used  in  the  specular  model.  In  addition,  the  depolarization  factors  DPO(B)  and 
DP90(B),  as  well  as  F(B)  and  G(BNP)  needed  in  the  volume  model,  are  defined  analytically. 
FUNCTION  FUNC  currently  yields 


DPO(B)  =  1.0 
DP90(B)  =  1.0 
F(B)  =  1.0 

G(BNP)  -  !  .0  r  BNP2  \ 

-  1 ;  - - 2  ) 

RCOSBNP  (COSBNP)2RPO  Qle'  e  2  SIGMA  i  ,  Q2  RHOV 

for  BNP  <  SIGMA 

/  BNP  \ 

-  (COSBNP)2RPO  [Qle  SIGMA  .  Q2  RHOV 

for  BNP  '  SIGMA 
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The  shadowing  and  obscuration  factor  applied  to  RCOSBNP  (measured  values  read  during  the 
input  phase  of  RHOPRIME  or  defined  analytically  in  FUNC)  is 


i  BNP  'TAU 
+  OMEGA  e 


OMEGA 


PHIEN  BEP 
OMEGA '  OMEGA 


(a)  SURFACE -PLUS- LAMBERTIAN  MODEL  CALCULATION 


RO  -  (N  l)2  i  K2  .  (V2  -  COSB)2  +  V3 
(N  -  l)2  +  K2  (V2  +  COSB)2  +  V3 

=  normalized  reflectance  lor  1  polarized  incidence 


where 


( V2COSB  -r  COS2B  -  ll2  -  V3COS2B 

- - o - cs - —  rlu 

(V2COSB  -  COS  B  t  1)  V3COS“B 
normalized  reflectance  for  •  polarized  incidence 


~j /  V4N2|- 


KnV_,  (N2  -  K2  -  1  COS2  B)2  -  (N2  -  K2  -  1  +  COS2  B) 

2 


If  H_  0  the  calculation  is  made  for  a  plane  polarized  source  (polarizai.on  angle  PSI).  The 
calculation  ignores  the  induced  elliptical  polarization  for  K  /  0. 


PSIELt  ATAN 


TAN  PSIOE  SIGN  tCOSATAN(Nl  -  CQSB)| 


polarization  angle  wun  respect  to  u,  £  reference  plane,  alter  reflection 

C  1  if  AP  0,  then  a  bidirectional  reflectance  Stokes  vector  is  computed 
C  AP  COSBKP  COSBDP  if  AP  '  0,  tlien  a  reflected  radiant  intensity  Stokes  vector  is 
computed 

nun  c  (K0  cos2pside  ♦  aau  sin^side 


j^COSBEP  COSBDP  - -  “ 

’  (RX1  COS2  PSI  PE  •  RX2  SIN2 


PSI  PEL 


UH2)  cj^^E^  cOSBDP(K°  COs2  P*,DE  K90  S,n2  PSIDF1COS2  (PSIED-WADE) 


111(3)  C  J-^lf^OSBDP  ,R°  COsi  PSIDE  ’  R90  SIN2  PSIDE1SIN2  (PSIED-WADE) 


VH 


i 

-i  j 

RCOSBNP  \  ,„r  1  1 

COSBEP  COSBDP  2  ^R°  R90)  +  2  (RX1  f  11X2 *  ' 


121(2)  =  C 

121(3)  =  C 

If  H  -  ±1  the  calculation  includes  the  phase  difference  and  ellipticity  induced  by  reflection 
for  K  *  0  and  is  an  exact  treatment  of  the  Fresnel  equations. 

SUBROUTINE  ELIPS1  (AA,  AB,  PSIDE,  H;  AA1,  AA2,  D)  defines  the  following  input  ellipti¬ 
cal  polarization  parameters:  the  amplitudes  perpendicular  (AA1)  and  parallel  (AA2)  to  the  D, 

E  plane  and  the  relative  phase  D  =  0||  -  01  of  the  amplitudes  of  the  major  (AA)  and  minor  (AB) 
axes  of  polarization  ellipse:  the  orientation  of  the  ellipse  with  respect  to  the  D,  E  plane; 

PSIDE;  and  the  handedness,  H. 

DR  —  0 1|  -  <t> 1  induced  by  reflection 

FOR  K  =  0,  DR  ---  0  if  COSB  <  COS  ARTAN(N) 

^  -77  if  COSB  >  COS  ARTAN(N) 

FOR  K  ,  0,  DR  =  *  ATAN  (-  2  Vy3<l-  COS2 BICOSB  - \ 

\  (1  -  COS  B)  -  COS  B(V2^  +  V3 )/ 

if  (  )  <  0 

DR  -  -  ATAN  (-  g.jVL!l^jg<g2B)C06B  - ^ 

\  (1  -  COS  B)  -  COS  B(V2  4  V3)  j 

if  (  )  >  0 

The  intensities  AIR  and  A2R,  and  the  relative  phase  of  the  parallel  and  perpendicular  compo¬ 
nents  of  the  reflected  radiance  induced  by  the  reflections,  are 

AIR  A1R0 

A2R  -  A2  •  R0 

DR  =■  DR  4  D 

SUBROUTINE  ELIPS2  f AA1,  AA 2,  DR,  AAR,  ABR,  PSIED,  HR)  defines  the  elliptically 
polarized  reflected  radiance  as  amplitudes  AAR  and  ABR  of  the  major  and  minor  axes,  the 
angle  of  the  ellipse  relative  to  the  D,  E  plane,  PSIED,  and  the  handedness,  HR.  PSIDE  - 
PSIED-WADE  is  the  angle  that  the  polarization  ellipse  of  the  reflected  radiance  makes  with 
the  normal  vector  to  the  OR,  D  plane. 

CHI  HR  ATANtABR  AAR)  is  the  parameter  used  to  d.fine  the  ellipticity  of  the  reflected 
radiance. 

C  1  if  AP  0 


§5§llloSBDpf<R0  -  R90>  C0S2  '-WADE> 

»  m 

eiHreiSBDP  i(R°  - R90)  SIN2  <-WADE> 


121(1)  -  C 
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C  AP  - COSBEP  COSBDP  it  AP  ,  U 
mm  r [AR  '  BR  WCOSBNP  _ 

111(1)  •  t  a  .  t>  pAcoPn  nAbiin 


111(4)  =  C 


"AR 

t 

BR 

A 

» 

B 

AR 

■+ 

BR 

L  A 

4~ 

B 

AR 

f 

BR 

.  A 

-+ 

B 

AR 

+ 

BR 

L  A 

+ 

B 

(RX1  COS2  (PSIPE)  +  RX2  SIN2  (PS1PE) 


-5  SIN(2CHI) 


RCOSBNP _ 1  /on  Dom  1,DV.  D  ' 

COSBEP  COSBDP  2  ^R°  +  R9°  +  2  ^RX1  +  11X2  ^ 


121(11  CLi^igloSBDP  I (R0  *  R90)  *  l(RX1  *  RX2! 

121(2'  ‘  C LSWeIwEDP  2‘,R0  R9W  C0S(-2WADE, 

121(31  C  jcWEWOSBDP  i  (R0  ‘  R90)  SINI-2WADE) 

121(4)  0 

(b)  VOLUME  MODEL  CALCULATION 

The  angular-dependent,  volume  reflectance  model,  Stokes  vector  is  given  by 
«'»<"  CDPM<PDP0>  Hp^IbDP  I00®2  PS1DE  DP90I1.DP0)  * 
Sl^PSIDE  •  DP90(1  - 01*90) J 

113(21  C  DMOil-UPO) ^P^BDPI00®*  W»*-OW<Kl-l*l>«)  . 
SIN^PSiDE  •  DPOO-DPSOllCOS2  2AD 


113(3)  C  ; 


2RHOV  E  G 


DP90(1  DPO)  COSBEP -COSBDP 


I  COS  PSIDE  •  DP90(1-DP90) 


SINT  PS  IDE  DP0(l-DPi)0)|SIN  2  AD 

1231,1  C  »W0ir.|lP»|  Sw'C^P  >P»U-DP0,  .  OPO.I.DPWII 
123(21  C  1)1*9011  DPO)  COSBEP -COSBDP  2  !DP9U  ~  ul>0l  COS1-2WADE) 


123(3)  C 


2RHOV  F  (J 


DP90(1  -DPO)  COSBEP  COSBDP  2 


'  IDP9U  -  DPO  I  SIN(-2WADE> 


a  here 


C  1  fur  AP  0 

C  AP  COSBEP  COSBDP  for  AP  *  o 


ml 


The  angle  of  polarization  of  the  reflected  radiance,  AED,  from  the  normal  vector  of  the  D,  E 
plane  is 


AED  =  ATAN 


DP0(1  -  DP90) 
DP90(1  -'DPOy 


TAN(PSIDE) '  SIGN(COt>  ATAN(N)  -  COSB) 


and  the  angle  of  polarization  referred  to  the  OR,  D  plane  is 


AD  -  AED  -  WADE 


SUBROUTINE  OUTPUT.  This  subroutine  prints  the  Stokes  vectors  for  the  bidirectional 
reflectance  (AP  -  0)  or  reflected  radiant  intensity  for  unit  incident  irradiance  (AP  *  0)  for  the 
surface  model,  the  volume  models,  and  for  the  combined  specular  and  volume  model.  Stokes 
vectors  are  printed  for  a  completely  polarized  source>for  a  completely  unpolarized  beam,  and 
also  for  a  partially  polarized  beam  (polarization  defined  by  the  input  parameter  P). 

In  addition,  several  calculations  are  made  with  the  Stokes  vectors.  For  a  bidirectional 
reflectance  (or  radiant  intensity)  Stokes  vector,  the  bidirectional  reflectance  (or  radiant  in¬ 
tensity)  for  a  receiver  polarized  1  or  II  to  the  OR,  D  plane  is 

A  +  B 

receiver  ^  -  — 2 — 


receiver  ^  ^ — 

where  the  Stokes  vector  is  of  the  form: 

The  angle  of  the  major  axis  of  the  reflected  radiance  and  the  percent  polarization  of  the 
reflected  radiance  are  also  given;  they  are 

at  -  *  ,TaX!  |C|  Qn  .  *  Qn  (looking  into  the  source,  AL  >  0  is  a  CCW  angle; 

AL-  i2-  AlANje!  -90  AL  90  AL  <  0  is  a  CW  angle) 


A 

B 

C 

D 


P 


x  kw; 


The  output  includes  TS,  PS,  TD,  PD,  P,  as  well  as  the  input  anu  output  values  of  A,  B,  PSI,  H 
from  the  surface  model  calculation  (if  the  input  H  =  0,  the  input  and  output  values  of  A,  B,  H 
default  to  1,  0,  and  0). 

SUBROUTINE  ELIPS1  (A,  B.  PSI,  H:  A).  A2,  DELTA).  The  basic  equations  which  relate 
two  specifications  cf  an  elllptically  polarized  beam  (A,  B,  PSI,  H)  and  (Al,  Bl,  DELTA)  are 


tan  or  Al  A2  0  i»  •  n  2 


tan  \  *B  A  for 


rt 

It 


-  77  4  •  V  ‘  77  4 


dl 


tan  2^  -  tan  2 a  cos  6 
sin  2  \  sin  2a  sin  o 


from  which  we  obtain 

,.:„2  n  sin2  2\  t  tan2  2^ 

1  +  tun 

or  equivalently 

cos  2a  -  cos  2  \  cos  2V' 

Subroutine  ELIPS1  determines  Al,  A2,  DELTA  from  A,  B,  PSI,  H 


LAMBDA  =  \  A2  +  B2 

If  B  -  0,  Al  -  LAMBDA  COS  PSI 
A2  =-  LAMBDA  SIN  PSI 
DELTA  -  U  when  0  r:  PSI  it  2 
DELTA  -  77  when  77/2  PSI  ‘  77 

Otherwise: 

CHI  H  •  ATAN(B  A) 

TI  COS  2CHI  COS  2 PSI: 

ALPHA  1  2  ARCOS(-TI)  if  .7  4  PSI  3tt  4 

=  1  2  ARCOS(TI)  if  PSI  <  77.  4  or  -  3t7,  4 
If  ALPHA  0,  Al  LAMBDA,  A2  0,  DELTA  -  0 
If  ALPHA  t;  4,  Al  A2  LAMBDA  V2,  DELTA  2CHI  if  PSI  -  ...-4, 

H -  2  CHI  if  PSI  -  3jj/ 4, 

If  ALPHA  r,  2,  Al  0,  A2  LAMBDA,  DELTA  0 
Otherwise 

TI  SIN  2CHI  SIN2ALPHA 
MU  ARSIN  TI 
Al  LAMBDA  COS  ALPHA 
A2  LAMBDA  SIN  ALPHA 
COSD  TAN  2 PSI  TAN  2 ALPHA 
If  COSD  0  DELTA  H  MU 
If  COSD  *  0  DELTA  H  (t7  -  MU) 

SUBROUTINE  ELIPS2  (Al,  A2,  DELTA:  A,  B,  PSI,  H).  Subroutine  ELIPS2  determines  A, 
B,  PSI,  H  from  Al,  A2,  DELTA. 


i>2 


LAMBDA  ]/a12  +  A22 

If  A1  =  0  or  A2  =  0,  then  A  =  LAMBDA,  B  -  0,  H  =  1,  and 
PSI  0  if  A2  =  0 
PSI  =  tt/2  if  Al  =-  0 

If  A1  =  A2,  then  CHI  =  1/2 1  DELTA  I,  A  =  LAMBDA  COS  CHI,  B  =  LAMBDA  SIN  CHI,  and 
H  =  1  if  DELTA  0 
=  -1  if  DELTA  <  0 
PSI  =  tt/4  if  CHI  <  tt/4 
=  377/4  if  CHI  >  77 / 4 

If  DELTA  =  ±77,  A  =  LAMBDA,  B  =  0,  H  =  1,  PSI  -it  -  ATAN  A2/A1 

If  DELTA  =  0,  A  =  LAMBDA,  B  =  0,  H  =  1,  PSI  =  ATAN  A2/A1 

If  DELTA  =  ±v/2,  H  =  +1  if  DELTA  >  0 

-1  if  DELTA  <  0 
If  Al  >  A2,  A  =  Al,  B  =  A2,  PSI  =  0 

If  Al  <  A2,  A  =  A2,  B  =  Al,  PSI  =  77/2 

Otherwise 

If  Al  >  A2,  ALPHA  -  ATAN  A2/A1 

CHI  -  1/2  ARSIN I SIN2 ALPHA  SINDELTA i 
LAMBDA  =  ITAN2ALPHA  COSDELTAI 
A  =  LAMBDA  COS  CHI 
B  =  LAMBDA  SIN  CHI 
H  --  ±  if  DELTA  l  0 

Part  1:  0  <  IDELTAI  <  tt/2;  PSI  =  12  ATANLAMDA 

Part  2:  n  2  <  I  DELTA!  <  n;  PSI  -  w  -  12  ATANLAMDA 

and 

If  Al  "  A2,  0  <"  DELTA  .  2  PSI  ..2-12  ATANLAMDA 

-  2  I  DELTA  <  r  PSI  -  ~  2  -  1  2  ATANLAMDA 

III. 2.  INPUT  DATA  FORMATS 

The  input  to  the  RHOPRIME  program  is  segmented  into  logical  blocks.  Each  block  is  ini¬ 
tiated  by  a  block  header  and  terminated  by  an  end  card  Blocks  may  be  in  any  order,  but  a  data 
block  is  assumed  to  precede  any  computation  request  blocks  or  scan  request  blocks.  If  a  block 
header  specifies  an  invalid  block  types,  all  input  up  to  and  including  the  next  end  card  is  ignored. 

DATA  TABLES  BLOCK.  The  data  tables  block  specifies  all  physical  characteristics  of  the 
materials  to  be  studied.  The  block  header  is  one  card  with  the  following  format: 


83 


Columns 

1-4 

5-19 

20-25 

26-80 


Description 

’TABL’ 

ignored 

maximum  material  index  to  be  expected 
ignored 


The  data  tables  block  is  itself  segmented  into  material  blocks  each  characterizing  one  ma¬ 
ternal  to  be  studied.  Each  material  block  is  initiated  by  a  material  header  and  terminated  by  an 
—  The  material  header  is  two  cards  with  the  following  format: 

Card  1 

Columns  Description 


Columns 

I- 4 
5-8 
9-10 

II- 20 
21-30 
31-40 

41-50 

51-60 

61-70 

71-80 


Columns 

1-10 

11-20 

21-30 

31-40 

41-5u 

51-80 


Card  2 


'MATR' 
ignored 
material  index 
n 
k 


SIGMA  (if  SIGMA  ■/  0,  RCOSBNP  is  computed! 
RPQ 

Description 

ignored 

« 

Q1 

Q2 

blank 


Following  a  material  header,  there  may  be  a  sei  of  p'  data.  If  present,  the  p' m  a  function 
of  f-  and  tiie  o . 's  must  be  in  ascending  order.  The  format  is 


Columns 

Description 

1-4 

blank  or  'ANGL 

5- 10 

ignored 

11-20 

"  N  (deg; 

n 

21-30 

Inn  u  n, 

31-80 

ignored 

WARNING: 

Each  material  block  must  be 

2 

|  COS  i' 


- - —  twwtvc 

block  must  also  be  terminated  by  an  end  card. 

CJ^MH^lION^EqUEST  BLOCK.  Ttie  computation  requests  block  contains  all  informa¬ 
tion  needed  to  perform  desired  computations.  The  block  header  is  one  card  with  the  following 


format : 


Columns 


in 


1-4 
5-19 
20-25 
26-80 

The  model  selector,  ISW,  is: 

1  —  if  specular  and  diffuse  models  are  desired 
3  —  if  volume  model  is  desired 
7 — if  combined  model  is  desired 

Following  the  block  header,  computation  requests  are  processed  sequentially  until  an  end 
card  is  encountered.  The  format  of  a  computation  request  is: 


COMP’ 

ignored 

model  selector 
ignored 


Columns 

Description 

1-4 

blank 

5-9 

ignored 

10-16 

source  zenith  (deg) 

17 

ignored 

18-24 

source  azimuth  (deg) 

25 

ignored 

26-32 

detector  zenith  (deg) 

33 

ignored 

34-40 

detector  azimuth  (deg) 

41 

ignored 

42-48 

polarization  major  axis  length 

49 

ignored 

50-56 

polarization  minor  axis  length 

57 

ignored 

58-64 

angle  of  source  polarization  (deg) 

65 

ignored 

66-72 

source  percent  polarization  +  100 

73 

ignored 

74-76 

handedness  of  polarization  (if  *  0,  elliptical  polarization  is  assumed) 

77 

ignored 

78-80 

material  index 

SCAN  REQUEST  BLOCK.  U  a  scan  m  the  detector  zenith  and,  or  azimuth  is  desired,  a 
scan  request  block  may  be  used.  The  olock  header  is  one  card  .  ith  the  following  format: 


Columns 


Description 


1-4 

5-19 

20-25 

26-80 


SCAN' 

ignored 

model  selector 
ignor«  d 


One  card  follows  the  block  header  giving  all  required  parameters.  The  format  of  this  card 


Columns 

Description 

1-b 

source  zenith  (deg) 

7-12 

source  azimuth  (deg) 

13-18 

initial  detector  zenith  (deg) 

19-24 

final  detector  zenith  (deg) 

25-30 

zenith  increment  (deg) 

31-36 

initial  detector  azimuth  (deg) 

37-42 

final  detector  azimuth  (deg) 

43-48 

azimuth  increment  (deg) 

49-54 

polarization  major  axis  length 

55-60 

polarization  minor  axis  length 

61-66 

angle  of  source  polarization  (deg) 

67-72 

source  percent  polarization  4-  100 

73-76 

handedness  of  polarization 

77-80 

material  index 

TITLE  SPECIFICATION  BLOCK.  A  title  may  be  printed  at  the  top  of  each  page  of  long 

form  output  using  the  title  specification  bl'c!:.  The  block  header  is  one  card  in  the  following 

format: 

Columns 

Description 

1-4 

'T1TL' 

5-19 

ignored 

20-25 

blank 

26-80 

ignored 

One  card  following  ' 

me  block  header  specifies  the  title.  The  format  of  this  card  is: 

Columns 

Description 

1-60 

title 

61-80 

ignored 

FACET  DEFINITION  BLOCK.  If  default  facet  definition  is  not  desired,  the  facet  may  be 

redefined  using  the  tacei 

definition  block.  The  block  header  is  one  card  in  the  following  for- 

mat: 

Columns 

De script  uni 

J  -4 

FACE 

5-19 

ignored 

2  U  -  z  5 

blank 

26  80 

ignored 

One  card  following  the  block  header  defines  the  facet.  The  format  of  this  card  is: 

Columns 

Description 

1-4 

blank 

5-9 

ignored 

10-16 

facet  area  (default  0) 

17 

ignored 

18-24 

facet  normal  -  x  (default  0) 

25 

ignored 

26-32 

facet  normal  -  y  (default  0) 

33 

ignored 

34-40 

facet  normal  -  z  (default  1) 

41-80 

blank 

8b 


END  BLOCK.  The  end  block  terminates  the  program.  The  format  of  the  block  header  is 
the  same  as  that  of  the  end  card. 

Columns  Description 

1-4  ’END' 

5-80  blank 

This  block  does  not  need  an  end  card. 
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Appendix  II 

INSTRUCTIONS  FOR  USE  OF  PROGRAM 
WITH  SAMPLE  COMPUTER  OUTPUT 

The  program  documentation  in  Appendix  I,  together  with  the  sample  computations  in¬ 
cluded  in  this  appendix  should  enable  the  user  to  (1)  modify  this  program  to  accommodate  the 
requirements  of  his  own  computer  and  (2)  verify  output  from  his  modified  program  by  com¬ 
parison  with  the  samples  given  herein. 

Note  that  the  input  parameter  values  shown  in  Table  VII  are  the  ones  with  which  the  pro¬ 
gram  has  been  run. 

Sample  outputs  presented  in  this  appendix  include 

(1)  a  listing  of  the  input  information  (Table  VII) 

(2)  the  computed  output  of  the  program  (long  form)  (Table  VHI) 

(3)  a  short  form  of  the  computed  output,  containing  only  that  information  necessary  to  feed 
into  a  computer  program  for  the  purpose  of  obtaining  plots  of  the  data  (Table  IX) 

The  three  tables  mentioned  above  appear  at  the  end  of  this  appendix.  All  of  the  sample  in¬ 
formation  is  keyed  and  labelled  so  that  elements  may  be  identified  easily.  However,  the  further 
descriptive  detail  below  may  be  helpful  in  studying  the  samples  given. 

RHOPRIME  Input  Listing 

The  following  items  appear  across  the  top  of  Table  VII.  One  line  2: 

n  -  real  part  of  index  of  refraction 
k  -  imaginary  part  of  index  of  refraction 
pxl  ,  cross  component  (2P|i)l  used  for  surface  mode, 

Px g  -  cross  component  (2pjj)J 
Py  =  volume  component  used  for  volume  model 
SIGMA  =  generating  function  parameter 
RPO  =  generating  function  parameter 

And  on  line  3: 

r  -  shadowing  and  obscuration  parameter 
SI  -  shadowing  and  obscuration  parameter 
Q1  =-  generating  function  parameter 
Q2  -  generating  function  parameter 

)o 

cos  6  tabulation  which,  in  this 
n 

case,  was  extracted  from  measured  data  and  determined  from  the  zero  bistatic  scan.  Alterna¬ 
tively,  such  a  tabulation  can  be  generated  by  use  of  agenerattngfunction  specified  in  the  SUBROU¬ 
TINE  FUNC. 
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Preceding  page  felank 


Note  in  tin'  sample  input  information  of  Table  VII  that  values  are  provided  for  p  . ,  and p  „ 

X  *  X  “ 

and  also  for  />y.  In  practice,  p  and  p^  will  be  used  or  py  will  be  used;  all  three  values 
will  never  be  nonzero  simultaneously. 

If  the  table  is  supplied  as  part  of  the  input,  the  parameters  SIGMA  and  RPO  are  set  to 
0  and  Q1Q2  1. 

TheP’(VA:  "a^a)cos2 

\  n  n  n  n  / 

the  computer  what  source-receiver  combinations  are  to  be  computed  and  what  model  is  to  be 


u  tabulation  is  followed  by  scan  request  information  telling 
n 


selected: 


rl 


r2 

"r3 

*r2 

*r3 

A 

B 

PSI 

P 

MI 

ISW 


c  for  source 

<p  for  source 

initial  »  for  receiver 

maximum  n  for  receiver 

size  of  angular  step  for  scan 

0  for  receiver 

0  for  receiver  (value  for  second  scan) 
size  of  angular  step  for  0 

semi-major  axis  of  polarization  ellipse  (normalized  to  1.0) 

semi-minor  axis  of  polarization  ellipse  (B  -  0  implies  linear  polarization) 

angle  of  source  polarization 

percent  polarization  (1.0  100 ,  ) 

material  index 


7  for  combined  model.  (When  volume  model  is  used,  set  p 


XI  PX2 


-  0.) 


Note  that  in  addition  to  these  input  parameters,  others  must  be  added  in  the  SUBROUTINE 
FUNG: 


DPO,  DP90  depolarizations  for  perpendicular  and  parallel  components  of  incident  beam 
t,  g  volume  model  parameters. 

For  the  materials  in  the  sample  listing,  values  for  DPO,  DP90,  f,  and  g  have  been  set  equal 

to  1.0. 


Cm nputer  Output  (Long  Form) 

As  exemplified  by  Table  VIII,  each  page  of  the  computed  output  corresponds  to  one  source- 
receiver  configuration.  Items  at  the  upper  left  are  self-explanatory.  However  it  should  be 
oorne  in  mind  that  MAJOR  refers  to  the  semi-major  elliptical  axis  (a),  which  is  taken  to  be 
1.0.  Since  MINOR,  which  refers  to  the  semi -minor  axis  (b),  is  0,  the  MAJOR-MINOR  combina¬ 
tion  implies  linear  polarization  with  polarization  angle  PSI  for  the  incident  beam.  HANDED  = 

0  whenever  the  polarization  is  linear  only. 


oo 


The  entries  in  the  three  main  columns  are  reilectances.  From  the  top,  the  first  four 
entries  in  each  column  are  the  surface  model  elements  of  the  Stokes  vector  which  describes 
the  polarization  state  of  the  beam  as  it  leaves  the  target: 

A  =  total  reflectance 

B  =  reflectance  with  receiver  polarization  angle  -  0  (perpendicular  polarization) 

C  =  reflectance  with  receiver  polarization  angle  -  45° 

D  -  reflectance  with  receiver  circularly  polarized 

The  second  four  entries,  still  in  the  surface  model  block,  are 

A+B 

— =  reflectance  recoided  from  receiver  with  analyzer  set  for  perpendicular  polarization 
A-B 

— 2~  reflectance  recorded  from  receiver  with  analyzer  set  for  parallel  polarization 

AL  -  angle  of  polarization  for  reflected  beam 
P  =  percent  polarization  of  reflected  beam 

Thus  far  the  first  two  blocks  of  four  entries  have  been  discussed.  The  foregoing,  as  pre¬ 
viously  stated,  apply  to  the  surface  plus  Lambertian  volume  model. 

The  third  and  fourth  blocks  apply  to  the  non -Lambertian  volume  model  and  are  to  be  inter¬ 
preted  in  exactly  the  same  manner  as  above. 

The  fifth  and  sixth  blocks  consist  of  the  sum  of  the  surface  *  volume  models  and  are  printed 
out  for  convenience. 

Note  that  in  the  volume  model  output  and  in  the  summed  output,  item  D  (circularly- 
polarized  component)  is  not  present. 

Computer  Output  (Short  Form) 

TTie  short  form  of  the  computer  output  consists  of  the  information  in  the  last  four  entries 
of  the  summed  output  (surface  <  volume*,  A^  ,  Ap^  ,  AL,  P  (see Table IX).  Moreover,  the 
data  are  compressed  so  that,  whereas  the  long  form  has  only  one  source -receiver  configuration 
per  page,  the  short  form  contains  a  complete  scan  in  one  block. 

One  scan  consists  of  four  item  numbers.  Preceding  each  of  the  first  two  item  numbers  in 
each  scan  are 

Wavelength  (1.06  pm) 

(0°) 

^(180°) 

*r  (0°  or  180°) 

A^B 

Tbe  first  item  number  in  each  scan  contains  — .  Each  output  entry  is  preceded  by  the 

scan  angle  —  i.e.,  0.0,  0.0266  means  that  the  reflectance  at  fc*  -  0  is  0.0268.  The  second  item 
A”  B  ^ 

number  contains  -y .  The  third  item  number  contains  the  polarization  angle,  AL,  at  each 
receiver  angle.  The  fourth  item  number  contains  the  percent  polarization. 

The  scans  are  in  the  same  overall  order  as  those  in  the  long  form  of  the  output 
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fcopni  3  APCS  TpOT  f  X  AP,  Y  1  1 
PHIF\  s  4.0 


IF ( (o'(»{ ' )-J°U ) ' 

.  L  T  .  ■  .  0  g  1 

.4*0.  Ao«(Yt?)-UP(?J ) 

.LT. 

||  .0001  ) 

ut Ty«N 

IF  (  A0<5(n(  1  J-JP(  1  )  1 

•  L  7  .  u  «  .)  p  0  1 

,A*u.  Aa^toCPl-uFlP ) 1 

.lt. 

J.0001  ) 

FtTO0'< 

IF  (  »h«(F  (  '  ).jt>(  1  J  1 

.lT.  n . j o g 1 

.  A  4.  u  .  A  oS  t  F  (  p  )  "jP  (  ?  )  ) 

•  LT. 

O.0O01  ) 

1  PkTo0N, 

jf(1  I  -  -«*!'  (aLuJ 
jo (P)  S  0.0 
0(*(3J  =  C7S°t° 

C 

Jl(1)  «  6  T  \  r  [jop  1  *f  u!?  t  F  uu„  ■  ) 

j'(?)  «  STNTtfOH)**!1 iFjPnl J 

O'  (  X)  *  C**SaU° 

c 

c * L »  PR''SX(r'c,u°.4i^i  ) 

P  *  VNO»*(W/1 ,N^1 J 


.04 


4*0 

4*1 

4** 

4*3 

4*4 

4*5 

4*6 

4*7 

4*8 

4*9 

400 

4«1 

40,; 

4«j 

404 

405 

4«6 

407 

40b 

409 

500 

5*1 

5**’ 

5*3 

5*4 

5*5 

5*6 

5*7 

5*8 

509 

510 

511 

512 

513 

514 

515 

518 

517 

518 

519 
5*0 
5*1 
5  ?i 
5*3 
5*4 
5*5 
5*0 
5*7 
5*8 
5*9 

530 

531 

532 

533 

534 
555 

536 

537 

538 

539 


C*u  ewnss(M^i  ,*c»N7) 

T 1  *  VN*RM(N2,N7) 

D*  *  *UT ( *  1 » N7 ) 

I  *  (  *N  .LT.  ft.*)°HTEN  a  1 . 5 7 0 70-ARCOS C -DN)~ 

H*TURW 

END 

SUBOOdTTN?  GFBNfC'1S*#CUSB*P»tns«EP,C08BNR,CnS*2.P*lPE»P31*E* 

1  4*DE#*P#A#AP#*»*R»M»HP»PSI01 

OTmFNSIIn  TABLE(S*0)»*TAB(500)#o(10)#9T(5) 
tnclC*L  MIDI  ,N0f)2,HnD'3#Mbn4 

E*U  I  VALENCE  (N00l,Mnu?,  Til#  (Mr»03,T31,  fHnD«»I4),  (  R(  1 1  ,N)  ,  («*(?)  #K  1  # 
1  C0(3),RV1), fR(4),P*P)# (O(5),RH0V). (9(6),HH0), 

*  (P(7),DP0), (R(B),0PO0), (R(9),FF), (R(l*),G) 

JNTFGPH  S4FTR 

WFal  121  (4),lll(4)tT23C3)»113f31#72«C3),ll4(3)#Tl?,T2?,N,K 
COMMON  MI  ,  I  4*.  «,  1  ABLE  .  i?l  .  1  i  i  #  1*3#  i  1 '3,  I?4, 1 14 
0*T*  MA8K/it*0*0*0*i/ 
feXTFRNAL  *I*N 

C *••••**•***•••***•***•***••*****•***••*«•*******•****•*•**** ********** 

c 

C  CwOnSF  model  AND  OETRTfcVE  D*  T  A 

C  “ . .  . 

C********************* a**************************** ******************** 

II  ■  l ANOf 4AS4, TSM) 

IX  z  I  AND f MA 3< , 8mF JO ( T SM , 1 ) 1 
14  *  1  ANDMAS* ,SmFTP(TSM,?)1 

cal',  getdat fR,cns«) 

O'!  10*0  1  *  1.3 

1*1  rn  «  *.* 

Ill'll  a  *.* 

1*3(1)  «  *,* 

113(1)  * 

1*4(1)  *  ".* 

1*00  114(1)  *  *.* 

c ••****•••*••**•••****. •*•*••*•••***•*•****••••*•*****•»***•••*•••*•*** 

c 

C  4rFC"L*H  and  0  I  F  Flj3t  MHuFC  ( M(JD  1  i  MHO?) 

c 

a********************************************************************** 

IFt.NHT.  MJDJ  .And.  .NOT.  Mnu?)RU  Tn  30*0 

c***«* 

c  LnMPoTt  FO£OnFL  C^EFF  T C TfcNTS 

c***** 

Vi  «  1'I7iD£(((D*1.01*(N*1.0)*rf*,^).t(''<“l.C)*(N”1.0)*F*4)l 

w  s  N*.N-**4-1.0*C1S«2 
p  a  sn4r('i,D*N«N***4*o*Oi 
V?  a  *«8T ( A"S( (°»n) *U.5) ) 
y3  a  A0ft( fP-w)*0.5J 

c 

k *  a  Vl*(f(w^-C'l5«J*(V2-C"»S*)*V3)/(fy?tCo3B)*(V?*Cj3B)*V})) 

c 

yfl  a  V2*C*1S0 

ROq  a  RD*((fv4*r34M?-l.n)**?*v3*C,T8«2)/((V4«C'38P2*1.0)**2* 

1  v3*CJ«B*)) 

1*1(4)  a  D , * 

111 (4)  a  *. * 

I*(M  ,NF,  0)61  Tj  2ft Oft 
AO  a  ft, ft 
bo  a  ft. ft 
081*  «  *. * 


10S 


auo 

hP  *  n.o 

buj 

C  *•  •< 

t  * 

c 

klAMfc  HTLAHT^ro  inyFL 

b<*4 

£****• 

T 1  a  AoS(bo^(°STOF)-l.b7o79b) 

3«b 

T  ?  =  Ad^(*bS(PSTDF)-9.71£lB9) 

but, 

I p £ t i  .r-T.  b.noi  .A*iu  T2  .r,T.  n.ooUPSTEft  *  Al*NfsoKT(r>ivir»E 

3U7 

1  (P9b,»0)  )*TANlPSint)*SlG.«fl,nSf  atan(M)  )-euf*b) ) 

butt 

I F  (  T  i  .  l  £  .  n.flol  ,U».  T  d  ,Lt.  b.bOUPSTfcn  ■  PSiuFaSIRNfCISf  ATAbt 

bU9 

1  (9)i-ru^d> 

b«!0 

APJ1  a  P$7£b.wAbt 

bM 

T  1  s  A 4 1 1 +  A°  )  1 

bS<? 

sT(i j  s  c">sm ) 

bM 

bT(?)  s  ST\fT1 ) 

bSa 

T  1  =  -WArtE-*AyF 

bbb 

y  1  S  (P0-KO0)*b.b 

bSb 

ST(tt)  S  (O0AkP0)*b.*; 

bS  7 

ST(U)  =  yi*rj<t(Tit 

b  bo 

b  T  (  S  )  a  V1»«1N(TH 

b^V 

c 

bbG 

1  F  ( 1  f  .cr.  b.bjRj  jn  ioiu 

bb  l 

v 1  =  »Hb/f C^bPtPbCU^Mbp) 

bb«! 

yu  «  t . b 

b*>4 

i;n  tj  lb^u 

b*>« 

c 

1  b  i  n 

vl  =  BmI.I? 

bbo 

yll  a  r  J«SdrP*C°SB  JP*  »H 

bb  / 

c  *  »  *  < 

r* 

bbo 

L 

£nMo0Tt  STjKtA  vtbinw 

b  b  9 

(.♦»*< 

b7v> 

C 

b?l 

t 

y^FnLAWT^Cj  bbyPCP 

b7d 

c 

>74 

1  b^b 

I?«!  a  io*i  ♦PA?J«vU*b.«; 

374 

ipirn  s  vi*5>Tr»)  +  i?^ 

37b 

r  V1*S7(U) 

b7t> 

I  P  1  f  4  1  a  V1*ST (b) 

37/ 

L 

b7  0 

C 

pn(.u,T/fj  btyPCF 

b  7  9 

c 

3*0 

IW  *  10*1  «r  j«(PSTMF  J**P*P»?*bJ»iPSTKF>**?J*VU 

3«1 

c 

3^ 

v  i  =  /  Z  "*5  ’  ►  a  i  3t )  *  •*  !♦  *1  *•'  C  »i  7  yF  )  *  »?»Psib  j  *  v  i 

3  M  -> 

I'lU'  =  t/ 1  ♦  A  » «r 

:,«4 

1 1 1  ( «  v  i  «jT i i ) 

b«b 

I »  1 F 4 3  X  yi*bTt ?) 

3  ®  0 

jb  Tj  4ao« 

3*7 

c  *  *  • < 

1  » 

3*0 

L 

tlL7PTlr*l  “Jbtl 

3h9 

C  *  *  *  i 

1  * 

3°i< 

,,T.  a  rJTKF.^,ialb«4 

b°l 

lr(OaTyF  .„T.  O.yWMbf  s  H«l">t*4.  l«l^9^ 

3<»^ 

lFlpbTyP  .WT.  v.ylHMbt  *  f-^nt*4. 1'U^VX 

3<>i 

A  A  a 

3°« 

A"  3  A. 7  I  t 3> 

b*>b 

C  *  C  ‘  FtTP,lfA**Ao.r«ibt,H,AAJ,M«<i,y1 

b«t. 

A1  «  A  A  1  •  A  A  1 

30/ 

A?  a  A  A  9 • A  A  P 

b«e 

t 

b<»9 

lF(AB*(rjAd-J  ,  |1  ,LT.  O.ObpyP  r  -i.mb9* 

10b 


600 

601 

60* 

603 

604 
60b 
606 
607 

605 

609 

610 
61  1 
61  2 
6 1  i 

614 

615 

616 
617 
616 
619 
6*0 
621 
62* 

623 

624 
62b 

62e 

627 

626 

629 

630 

631 

632 

633 

634 
63b 

636 

637 
e36 

639 

640 

641 
64* 

643 

644 
6«b 
6«0 
64  / 
640 
644 
650 
6^1 
65* 

653 

654 
65b 

656 

657 
656 
659 


IF(*b5(P(J3B)  ,LT,  0.001HJ9  s  n.o 

c 

1P(K  .NF.  0.01GH  TO  *640 

bPTAb  =  C"»Sf  ATAN(N)  )  . 

IF (COSp-SFTAb 1*0 10, 2020.2030 

*010  QP  «0.0 

GD  TO  *«60 

20*0  JP  a  -1.570796 
Gn  To  *06« 

*030  1)0  a  -3. 141593  - 

Gr>  TO  2060 
C 

*0')0  t3  •  2.o*4q07(v^)*(  i , 0-CQ9B2) *Cf>36 

T4  *  (1.0-Cr'S521**^-cnS52*(V**V27V3) 
1P(485(T4)  .gt.  O.OOOODGn  TO  2050 
0»  *  -1.570794  ... 

an  to  *060 

c 

*ObO  j 1  «  T3/T4 

lF(Tl  , LT.  O.fljOH  a  -3.l4lb43+ATAM(-'t ) 

1 5  {  T  1  .  G£  .  0.0)09  a  -A T *N ( T 1 ) 

C  ■ . . - 

*060  ()P  a  flRTQ 

AIR  a  A 1 *60 
A 2r  s  A2*B90 
a  A i P  s  Ago) (AIR) 

AA2P  a  Agar  f a2r ) 

CALL  FLTPA*f AA1P, iA>R,gP, A A» , AoP , PS T£0 , MR ) 

A°  a  AAP*tAP/r A*B) 

6®  a  A6P*A6P/ (A*d1 

IF(PSTto  ,  GF  .  3. 14l593)PSItO  a  PSTfcO-3. 1 4 1 593 
p«10  a  »sTEO-«AOt 

IF  (®ST  J  .OE.  3.l4lb°3l1»3in  a  P6TD-3. 14j593 
Ip(pSTO  .  i.  £  .  0,0)05Ti)  a  PSI0t3. 1  41  593 
T 1  a  PSTg+PAin 
T2  a  I.970 7*6 

I5(»b5(AAP)  .r,T.  O.O0ODT2  a  HR»AT*N(  AtJP/AAP) 

T 2  a  T *♦ T 2 

C 

1 5 ( A  p  .OT.  0.0)OJ  Tn  ?070 
vl  a  Pun/ rc''S5k2*PU5bop) 

V2  a  1,0 

an  T L)  *060 
C 

*070  v 1  a  Suit  A p 

V  2  3  r  jAdFpaC'TS"  JO*Ap 

c 

c  cnMp<.  Tt  3T04tA  verrnp 

£>*«•* 

C 

C  a 4 t 2 f )  b'lgPcP 

C 

2«60  l?2  a  f»viTR**)av**0.5 

I2lfll  a  V1«(90»R601*O,5tT*2 

121(21  8  Wt*(P(;-R90)*0.4*Cu.3(-ai05-PAOt) 

1  2 1  ( 3 1  a  Vl*(P0-9901*0,9a3JN(.aAt>F«4Ar>tT 
C  PO'L»RT*fo  S1U6C5 
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Oh  1 
662 
663 

669 

66b 
666 
66  7 
66B 

669 

670 

671 

672 

673 
679 
67b 

676 

677 

676 

679 

0*0 

6*1 

6*2 

6*3 

6*9 

6*b 
6*6 
6*7 
6*6 
6*9 
6<*0 
6**  1 
6  02 
6°3 

6°9 
6°b 
6<»6 
6«7 
60« 
609 
7«y 
7  0  1 
7«2 
703 
7  0<4 
70b 

/o0 

707 

/oe 

7  «9 
/U 
71  l 

712 

713 
719 
71b 
716 
71  7 
716 
719 


C 

112  a  M91*CnSfP,Ipfc1»*2*K#2*STNfH0IPfc)**2)*V2 

C 

va  s  vi*r*o+"«i 

mm  =  v3»n,i 

mm  *  v3*cnSfTl)*cnsfT?) 

111(31  S  V3*STN(T  n*Cns(T?3 
Ill(a1  a  V3*STN(T?} 

C******************************** ******  Oft************************ 

c 

C  VHLHHF  MJOEI,  (MO03) 

c 

C** ********** ••••*****••*••** ******** ******•••**••* •****«•***•**• 

3000  IF(.NnT.  «J031&'1  TO  90 0 0 

VI  =  1.0 

IF(AP  .  GT .  O.OJV1  s  tflSOtp*r0SBOP*iP 

Vi  a  ?,0*vi«WM(jv*Ff.*(i/(CU9BeP  +  CnsP0P) 

*F0  =  O.y 

T 1  a  *B«(Ao*(PSIUF)-l.b7y796) 

T?  s  *BS(*aS(PST0F)-9. 7123BO) 

1 F  ( T 1  ,r,E.  o.oci  t  2  .Rt.  o.O01)«l)  T"  3010 

AFd  s  H«UO£ 

Gn  TO  3030 
C 

3010  IF(r>K9u  .it.  0 .  no  1  .AMO.  npo  ,LT.  O.909)(,n  TO  3020 
*Fu  a  1 . b 7u 79 
On  Tj  3030 

30*0  v?  s  Ouotopo  •  (i  .o-uPvoj /(OHou,(  t  .n-npnj ) )*r»N(pSifitl*nir:N 
1  (c'’S(AT*»i(9)')-ruobi 

AFd  s  ATA9(V21 

c 

3630  *0  S  AtO-wOt 

Cl  a  rj9(OSTDF) 
ci  a  ri«C1*opou 
«,1  a  <UV(0STDF) 

31  a  M*S1  *opn 

Vi  *  Vl/(fl.0*i)Pu1*Op901 

(••♦** 

c  CnMPuTt  S TD^t * 

C***** 

C 

c  pn^ArfT^PO  SOyOt^ 

C 

113(11  a  VI  *(r  j  ♦  (  1  .  n*npn  «f  i  .i;*uuvn)  ) 

113(21  *  vl*(Fl  •<  1  .  O.rVOjTSt  «(»  .  v-i>P90)  )*ry9(»y*An) 

113(31  =  vi*(f  i  *(  1  .  o.nyn)  as’ *f  1 . 0-GP90)  AyAAflj 

C 

c  u^PClaxT afd  b^yPCF 

C 

.1  *  V  1  *  0  .  b 

1*3(11  a  VI*  (OyOy  *(  1  1  .  n*rv6y)  ) 

1*3(21  a  Vl*(lHOy-oPy1*roS(«.,*ilF-'*AOt  1 
1*3(31  a  vi * (oy»o*yPo 1 •*! 'u -*AyF.*AOt ^ 

{.*•*«••«• a**************************** ..*•••*..••••••• ••••*•••••» 

C 

C  T  i^Ti  OF  *  Cc  "JOitl  (Mncl9J 

r 

w 

C •**•*•*. *•••••*••*••«*•.•••••*•.••••«•*•«••*••••••*•••**•**••••• 

9  0  0  0  1F(.VH.  9  J6u  T  MF  T  '19  *1 

kFjUrfM 
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720 
7?1 
7  22 
723 
729 
725 
7?6 
727 
72b 

729 

730 
V  3 1 

732 

733 
739 

735 

736 

737 
73b 
739 
79u 
7«1 

792 

793 
799 

795 

796 

797 
79b 
799 
7«;o 
7M 

752 

753 
7544 
755 
75b 
757 
75b 

759 

760 
7M 
7*2 
7*3 
7*9 
7*5 
7  *o 
7*7 
7*8 
7*9 
/T  ; 

7  7  l 
(?2 
77  3 
7  744 
775 
770 
777 
77# 
779 


£ND 

SIlbPJIIT  TnF  OUTPUT  (TS,PS,TO,PD.P,AA, AR,bB,RR,PSJ.P5in,M,MR,LABFL» 

1  AP) 

DTmFNSI'IN  Ar3),»(3),Cf 3),  »E.5(  3,  3),  AFpji,  3),  ALt3.il/2Pf  3,3V 
DIMFNSir»N  l*BFL  f  I5),0f  1),  TABLF(500) 

L^fiTCAL  MIDI  ,9ur»i,  (-nD3,9Q(l9 

RFAl  T21 (a), 11if9l,T23(3), 113731, 729 (5), 119(3) 

1NTFI.FN  S9FTW 

cnM9jN  E.  121,  IJ  1,1 23, 1 13, 129. 119 

eruivalfnce  (mooi.moop, ti V , t'nr)i»7,T3), tmoia, t«i 

data  M A*K /zn 0 00 <109  1/,?tPU/0.0/, (INF/ 1 .0/ 

c  ***** 

C  Formats 

C***** 

100  FORMAT ( T 1, 1  ft* , 'PEFlFCTANCF ' ,2*,9L1 ,2*. 15A9//21X, 'ThFTA ' ,6*, 'PHI • , 

1  13X, 'PJLAP^t"1  . 13*,  •UNPuLAPITEr''',  12*, 'PAOTTAL  Pui-.'J 

1 1 0  FORMAT ( 7 1 , 1  ft* , < 7NTENS7  TV  , 2*, 9L 1 . 2* , 1 5A9//21 * , • TMFT A • ,6*, 'PHI • , 

1  13X,'O3LASl7E0',l3)i,(ilNP0LA(>I?fcfl',l2*,'PARTTAI  Pul.') 

l2o  FHrmaTII 1 v, ,SnUsCF,»2f9X,F6.2),7(10*,Fi3,5)/11X,,D£TteTrtR,,7X» 

1  FA.P,0X.FA.2,3( 10X.Fl3.5i /l ix, 'PEBCFnT  PQLAP. * ,*X,F6.2,3tlOX, 

2  E13.5V37*,3(10V, Fil. 5)) 

130  F  nWM A  T ( ?3  V , 'lN',8*,,U"T,/11*,,MAjn(«(,1x,279X,P6.2),t(!0X,Fl3.5)/ 

1  11X, 'MiMjR',iv,?(9x.F*.2J,3tlft*,t1i.5)/1llt, 'P*I '.3*,27o*,F6.21 

2  ,3(107, P13.5)/11*. 'MANDFu*,P(a*,FS.ft 1,27,3(1 07, F13.S)//) 

190  FORMAT( 7 (377,3(1 07,F13.5)/)/9f 37*,3f l « X , t t 3 . 5 1 / 1 // ) 

PA  =  181(9) 

IT  =  1 

11=1  ANl)f  HAS*,  7{H) 

IT  =  I  A  7  HA  S*  ,  5lFTP  (  7  S*,  1  )  1 

la  X  |  A*'D  7  HA  s*  »  ShF  TO  (  Ts*,  2  )  ) 

IF  ( .Nil .  -JH  .ANft.  ,  N  u  T .  wnD?)r.u  TO  Po^O 

A  7 1 1  X  “A*I11fl1*7l.0-P*)*121711 

8  7  1  1  =  PA*I11721*71.U-PA)*1P1721 

C  F  1  1  =  OA*1117il*7l.(i-pA)*l?l7il 

L)  F  1  )  8  PA*H1791*71.0-PA)*121791 

AFS71.11  x  F  J,  «  1711*711  (?))/?.  ft 

aFS72.11  s  tI2ltli*T2t(2)l/2.ft 

AFS(3.n  =  7A711*»(1  Jl/2.0 

AFp7i.l1  x  711 1711-711 (2)1/2. ft 

AFp72.ll  *  '121(11-121 (2)1/2. ft 

AFP73.H  x  7A71)-P(1))/2.0 

Al  (  1  ,  1  )  X  n.OV99«t'> 

Ir  (  7  1  1  (?)  ,*t.  9.(1  ,u».  m'3)  .NF.  u .  0  1 
1  AL7l.ll  X  A1AN?(T11  (T),iH721)*2*.*9"b 
At  (2, 1 J  x  ft,«9*>9PtP 

Ir  (  721  (  *  )  O.ft  .UP.  I  ?1  7i  1  ,*F.  U.ol 

1  *Lf2.11  *  *79*2(72)  (3),l?i72l)*2*.*9*b 
»'  (3,1)  x  n.aR9voto 

IF(P(7)  .‘’t.  9.  (1  t  (  1  1  v.  .01 

i  * l  1 3, 1 1  *  a  :  *  .P(r  i » ,  ,.m*  i ) )  *2P.*9«b 

pp(l,»>  X  '  .9  (  r  *  1  4  (21  •  T  1  »  (  2,*4  1  I  f  >1*71 1lT)*  1  1  1  7  44  1  *  7  n  (9)  1/ 

1  711i1)*ifto.(, 

pP  (  2  ,  1  j  s  <j»l71Plf2)*72U2l**?l731*T21(3)*l2l79>*721  (9)1/ 

1  721 (1)*lft0.O 

PP  (  3 ,  1  )  «  Hy»i  7  o  7  1  1*B(  1  )♦(.  71  l*r  ( 1  )*y  7 1  )*(1(  1  )  1/A(  1  )*l«o  .0 

2ft Oft  1F(.mO|.  t.j  >ftoft 

*721  •  BA*113711*71.0-PA)*1P3(11 

b(2 )  x  »A*Il3f2'T7l.0-P*)*12372l 

C(pl  x  09*113(31*7 1.0-P*)*123(31 

AFsfl.21  *  (li3(11*Tl3<2)1/2,ft 
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7*0 
7* 1 
7*2 
7*4 
7*4 
7«S 
7*6 
7*7 
7*B 
7*9 
7<»o 
7°1 
7p2 
794 
7  9a 
79b 
i"B 
7«7 
7*6 
709 

e*o 
6*1 
6*2 
6*4 
6*9 
6*S 
6*6 
6*7 
6*6 
6*9 
610 
61  1 
61* 
6  1  4 
614 
61b 
6  1  6 
617 
4  1  6 

M  9 
6P0 
0?l 
6?* 
o?4 
bp« 

e?o 
0  J>7 
e’e 
o3* 

B  Vj 

*M 
*3* 
6*3 
6  Tm 
e3b 
fl  3b 
83/ 

6 '6 
634 


*Fs(2,2)  s  fl?3m*T*1£?»/?.n 
AFS(3,2)  =  (  A(2)T*(2)  )/?.* 

APpi’i^i  s  fii3m-u3(?n/?.o 
A*p(*,*i  3  fI?4fn-T*M?))/?.n 

AFp(4,,»)  »  (A(2)-«(2)1/2.0 
A I  (1,?)  s  *,99999t9 

IF(Tl3(2)  •  klE  .  *.*  .Up.  113(3)  .NF,  o.O) 

1  a  AT*n2(T13(3), 113(2)1*2*. 64*6 

Al (2,2)  s  * . 99«99£9 

1F(T23<2)  .ME.  O.n  .OP.  123(3)  .N>.  0.0) 

1  «L (2,2)  s  ATAN2(T23(3),12i(2))*2*.64*8 
AH?,?)  s  n.09999£9 
I  F  l  *  (  2)  .ME.  O.fl  .  1)9 .  cr*)  .NF.  0.0) 

1  Ai.73,*)  3  ATAN?(r(?)  ,6(2)  )*2».64*B 

Pp(l»?)=S*HT(Tl3(?)*I13f2)»Tl3('X)*I14(i)^/113(l)*)0*.* 
pP(?,?)sS*4T(T<;3(?j*i?3f2)+T*3(3)*l?3F3))/l?3f 1)*10*.0 
pp(3,2)sS04T(*(2)*B(2)*C(2)*C(2))/A(2)*l00.0 
3*0*  I f  1 . n*>t .  *<3*4)00  TU  a*u* 

Hafl)  s  Til  (1 J  +  I130) 

124(1)  3  T21 ( 1 ) ♦! ?3 ( 1 ) 

114(2)  3  Til (9)*I13(*1 

124(2)  s  T21 (9)AI?4(2) 

114(3)  s  TU  (3)*I13(3) 

124(31  3  Tgi (3)^1»3(3) 

A(3)  3  PA*114(11+(1.0-PA)*1?4(1) 

6(4)  S  PA*Il4(*)+(1.0~pA)*124(£) 

C(i)  s  PA*I 1 4(31+ ( 1 .O-PA )*i?4(3) 

AFS(»*3)S  t  T 1  a  C 1  J*I14(*) J/^.O 
AFSf*.b)»  (T*4(1 J+i?4(*) ) /*.0 
AFS(3.3>  3  (A(4)**(3l 1/2.0 
AFH(1.4)  =  (Il4(l)-T14(?))/?.* 

AFP(<|,41  3  (l?a(H-T24(?))/?.n 
AFP(4,3)  =  (A(4)-*l3))/2,0 

Al  (  1  .  3  )  3  0.94999*9 

1F(T1«<>)  ,Mt.  *.*  ,up.  Jta(4)  .NF.  U . 0 ) 

1  »L(1.3)  *  »TAN2(TlU(X).11a(*))***.64*6 
Al <2,31  3  0.9999939 

1F(T*«(2)  .ME.  *.*  .(jO.  1 2a  (  4)  .NF.  o.u) 
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Al(3,3)  3  0.9949939 

I  F  (  *  (  3  1  ,9£.  0.0  ,jQ,  ufi)  .NF.  0.0) 

'  ‘L(3,i)3  AT**',;  f  0  (4)  ,R( \)  )  *2b.0'IMA 

pP(’#'J3S*4TjT,a(2}*iiaf(iitTiari)*ii4(4))/li4(i)Alon.* 
ppt2,3j3a''4T(T24(»j*12a(*)^T2U(T}*I24(i))/l2«(n*lo*.* 
PP(3#3)rb,)4T(a(Tj*oT4)Tr(3)*c(4)l/A(3)Alu*.* 

C  A  Cl  AyTI')(LA*eL,Tu),«Fbfl,i),Ato(i#'»),AI  (1,3),pP(1,3),2> 

4*0*  13(AH  .it.  *.  *)  MKT  1  3  (^,  1  o*  J  T1  ,,-irtoi  ,MUOi,Mnj3,Mor>a,E*pFE 

IF(*P  •  *  1  .  *.*)woT  13  (  A,  1  1  0  J  1 1  ,Mnot  ,MUrti,Mnl)3,**Ur'4,CAoFL 

I  F  (  ,  N<)  I  ,  «J*t  .an*.  ,*'OT.  mHu2)OJ  J  O 

opITt(i.l>0)T*,pb,  111(1),  7*1  l  1  ),Afl),Tu,pl),Tl  1  (2l,  I>1  (*),  "  l  1  )  ,  PA, 

'  mo), i*i(3), r(i),ii:(4),T«>im,j(i) 

oplTE(3.no)A».,A4.  (Afc*  (T,1  ),1*1 ,4>,«o,rtP,(  AFP  (1,1),  T*  1,1),  K41, 

1  p3T0,(*Lrl.l),T«1,3),H,H9,(PP (1,1 1,131,3) 

lF(.„'lf.  -4*3)4^  TJ  4*1* 

*°lTt  (3,  100)113(1),  T21(1),A(2),Tll(P),X?4f^),(»(>),11>(3),T(ein), 

’  F(»).(At*lT,>),I»l,3),(AFpfl,2),T*1,3),(ALfl,2),T*1,ll, 

?  f*9(T, >),!*), 3) 

4*1*  1 F  (  . ;,n f  .  M3>)4)h31i|4*( 

*°ITt(3,  1*0)114(1),  T24U),*f4),Ti<i(?i,i?af^),*liJ,xi«f3),T2«Cl), 

1  r(').(*t*(T,^J,lBl,4).(AFPfi,3),T»1,l),(*Lfi,4),T»1,l), 
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8«! 

842 

e«3 

B44 
B4b 
846 
BA  7 
8«a 

849 

B^o 

BSl 

6S2 

8*3 

8S4 

8Sb 

bSb 

8S7 

B'iB 

B^V 

a*o 

8M 

8*2 

8A3 

8*4 

8*>b 

bb6 

8*7 

8*8 

6*9 

H’O 

B71 

672 

6’3 

87« 

6?b 

B7fc 

677 

rt7« 

079 

8*0 

8*1 

8*2 

8*3 

8*4 

8*b 

0*8 

8*7 

8*8 

8*9 

6°0 

8<»1 

B°2 

d«3 

6«4 

6«b 

8*8 

8*7 

8«8 

8«9 


?  fPP<T,3),I=l,31 

WFT'lkN 

4*2*  1*1.  NOT,  «JP3)GP  TU  4*30 

«9ITtf3,l*01T«i,PS.I13fn,T23(1),Af2)»TD,Hn/T13(?),l?3f27,R(?).P*, 
1  T17(7), 1*3(3), C(?) 

^oiTEfJ,  170)CjNE,ZF*P,  f  AFSU,2).Te1,3),ZFkn,7fcBi),  ( *EP C T , ?) . I  =  1 . 3 1  , 

1  PST,7E0j,C*ufi,2), Tsl , 3) ,2FRn,7t9Q. C PP (I , 2 1 ,  Ta  1 , 3 ) 

IPC.NHT.  MJ641WP1  '  IW81 

8»lT£f3.140Ut4fll#?24C1}.*f3),Tl«(?J.I?4f2),«(3j,ll4f37,T2a(3), 

»  C(X).  {*t1(T,3).I=1.3),  fAFPf'l.3i#Tsl  ,3),{*L(  I.3l,T»1 ,3), 

?  fHP(T,T).I»l,31 

M  *  T 1 1 K  *i 

4030  IPI.NHT.  «3r)4)KFTMHN 

*PlTfc T3, 1 *0) T *,PS, 11 4(11,124(1), A f3),T0,pn,T 14 (?),l?4 (21,* (3),PA, 
1  114(3). I?4(3),C(3) 

8»lTEf3, 130iUMt.2*kn, (AFS(I.3j,T=i,3),2FHn,7£Oij,(4tP(T,3),I*l,3), 

1P*I.ZFHp,  (Al(T,3),I  =  l,31,7fc»i),ZFRn, (PP(T,3)»i=l»3) 
r(F  T'lHN 
£►'0 

SMbPU'ITTnC  AuvmtTlTLP.  *1  .  Yl  ,  V2 , Y 3 , V4 , I CuDE 1 

CnMHU»!  /CM*T/0S1  ,  oui  ,«fc.TA,bFl  A6,kPU,CnslNiF,Sir,M,PHTfc7',«tP,TS 
OTmFNMPN  TTTI  fcC31,V(9,48l 
OAT*  9t/1.08/ 
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c 

C  SI 'bp J11  T  tNc  rjMToji 

0 
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c 
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c 
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4«u 
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9M 
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c  *  *  * 

Si"* 

t 

904 4 

c 

1 te  jut 

90b 

c 
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c  *  *  * 
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3000 

cnsTiMjr 

90fl 

c  ***** 
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c 

HFFLtr f ANr£9 

91  0 

c  *  *  *  *  * 

91  1 

po  =  Pdl *37 . 2°b7 / 

91-2 

P3  a  P3i*S7.£9b77 

913 

O')  3010  1*2,3 

919 

I  *  =  TS+l 

91b 

*PIT£f«,  lions,  TP,  wl  ,TS»PS,Pj 

916 

*OITE(4,1?01C7(I , Jj.Xfl.dl.Jsl.TPT 

91  7 

3010 

C'lNtlMUF 

916 

c  *  *** 

* 

919 

l 

A 0- li «  fc  UP  OjLAOlJATjnN 

4?u 

c* **** 

9?1 

13  a  T  b* 1 

9?* 

*PITfcfn, 110>I9, TP 

9?3 

oPlTfcfa.noMVl  J  ,  H.XfU.J),  r=1,TPT 

9?4 

c  ***  * 

• 

4?b 

c 

pppot^r  P"*L*t<T/ATTOW 

9?6 

c  *  **• 

* 

9?/ 

13  =  T  S* 1 

9  ?6 

*°lTtf«4, 110)19, TP 

9?9 

*PlTfcf4,luo)(»(1,  1) ,*fb, J),  la  1 , T P 1 

930 

JP  a  r\ 

431 

KPTIIpK! 

432 

t*u 

933 

i>ll0PjllTT<i,F  FLTP3lfA,B,P31,rt,Al,A2.uFLT*) 

43« 

PP  Al  1  PH 

93b 

lp(pi>T  .Lt.  0.0  ,*lp.  Pat  .t,7.  3.iaiS43)C,j  JO  n0no 

93e 

t_  *  M«lb  A  3  3  ,B]  f  A*A*d*&) 

37 

c  *  *  *  * 

♦ 

*  <tt 

c 

jPoFnfhaip  r a 

939 

c  *  *  *  « 

* 

4  <J  0 

1F(“  .t»7.  9.0)001  TJ  1O0O 

1 

A’  a  Aa9(r3a(o<,Tji*i  A“r>OA 

94<; 

A5  a  Jb1(H‘!(P!,T))«|  A**Dn* 

4«i 

pfl  T  *  «  0  .  U 

*'*« 

lr(pbT  ,„T,  1.67j79A)Otl  j  A  8  3.1u1b93 

9ub 

PP  THh*' 

4  9e 

i  *  *  *  * 

• 

4'i  7 

c 

t'  lTptkai  r  a  a  t 

9/J* 

* 

949 

10U0 

C  u  I  a  h*ATa«(W/A) 

9*0 

I<  s  Arta(Pj«t?.n.r„»)*(,ni>f<.u*p3|T) 

**  1 

lc(PaT  ,„T.  0,/»jl4«  .A*u.  p^T  .1,7.  i .  33nl  49  j  iL»HA  8  .  b*  APcri*  (  .7  1 ) 
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IF(°i*T  .tF.  0.7*3'4O  .  j  p  ,  °6  T  ,(,P.  *.39019*)ALP„A  »  .bAAPCOi-Sf  7jt 

4^3 

iF  l  ItPol  .»P.  0.0)4')  Tj  1  0  J  0 

99(4 

41  a  1  A*'*'** 

4«.b 

*P  a  0 . n 

4C6 

iiP(.ta  a  0.0 
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9«6 

101« 

lP{*b9(Aw3pA-0,703A64)  ,j7.  u.uOoljPH  t 1  1W?0 
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Al  a  |  A«tOA/l.«tl9?l'l 
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960 

961 

962 

963 
969 
965 

96fc 

967 

968 

969 

970 

971 

972 

973 
979 

975 

976 

977 

978 

979 
9*0 
9*1 
9*2 
9*3 
9*9 
9*5 
9*6 
9*7 
9*8 
9*9 
9*0 

901 

902 

903 
9044 

905 

9°6 

90/ 

90t» 


c 

1*2* 


C 

1*3* 


A?  a  A 1 

liFLTA  s  2.0*CMI 

OFLTA  -  M*5.101S91 

K  r  I  UK  N 

lP(*i<*(AL*HA-1. 5/0706)  .ST.  0.0001  ) Gu  T<7  1030 
A 1  s  0.0 

A?  *  l  AM80 A 

U*tTA  =  0.0 

9*  TURN  —  -  - 


Tt  =  A85(SIM(?.fl*CHT)/srNC2.0*AI P«A ) ) 

IP(Tl  .r.i,  1  .  o  j  r  i  a  i.o 
M"  a  AH5I9(T11 
Alai  AM8Da*C'ISUI>HA) 

A?  a  t A«8*A*STN7*LPHA) . 

cr'sf'  a  7AM(?,0*P57)  /  TAN(2.0*ALPHA) 
ip(ro*o  .nr.  *.*)*tt.T*  =  h*mu 
IP  (7*0*0  .  L  £  .  *.*)*£*  T  *  a  H*(3.  l«1593-Mii) 

MPTIIMH 

C 

C***o* 

8*0*  oPITt (3, 1 *0)P91 

ST)P*TC///I  *****°ST  nUT  r»F  HAN5E  -»',M0.3) 
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C 

C* 
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» 
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8  a  0.0 
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1* 
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c 
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